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Abstract. We survey some results on scalar curvature and properties of so- 
lutions to the Einstein constraint equations. Topics include an extended dis- 
cussion of asymptotically flat solutions to the constraint equations, including 
recent results on the geometry of the center of mass of such solutions. We also 
review methods to construct solutions to the constraint equations, including 
the conformal method, as well as gluing techniques for which it is important 
to understand both conformal and non-conformal deformations of the scalar 
curvature. 



1. Introduction 

The mathematical influence of Richard M. Schoen can be measured in many 
ways. His research has fundamentally shaped geometric analysis, and his results 
form many cornerstones within geometry, partial differential equations and general 
relativity. Evidence of his influence includes the large number of his students who 
continue to work in these areas. As two of these students, the authors of this 
contribution are exceedingly grateful for Rick's mathematical insight and generosity. 
Many of the results discussed in this article owe a great deal, directly and indirectly, 
to his continued mathematical vision and reach. 

By the second half of the twentieth century, the theory of nonlinear partial 
differential equations (PDE) had developed sufficiently to allow more sophisticated 
analysis to be brought to bear on geometric problems. One such early triumph 
was achieved by Choquet-Bruhat in 1952 [53 , when she exploited the underlying 
hyperbolic nature of Einstein's field equation to establish the local well-posedness 
of the Cauchy problem in general relativity. This breakthrough ushered in an era of 
intense study of both the (essentially hyperbolic) evolution problem for Einstein's 
equation, and of the set of allowable initial data for this problem, which must itself 
satisfy a nonlinear elliptic PDE system, the Einstein constraint equations. This era 
continues to the present day, with many remarkable recent advances and continued 
interest in the field from both physicists and mathematicians. 

Another watershed for the mathematical development of relativity is the cele- 
brated work of Rick Schoen and S.-T. Yau from the late seventies on the Positive 
Mass Theorem. Not only did their work employ serious tools of geometric analysis, 
including partial differential equations and geometric measure theory, to resolve 
a question motivated by gravitational physics, but they also established a link 
between the positivity of the mass of an isolated gravitational system and the re- 
lationship between positive scalar curvature and topology, a topic of interest to a 
broad range of mathematicians. In the early eighties, Schoen brought the Positive 
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Mass Theorem to bear on the resolution of the famous Yamabe problem, provid- 
ing more evidence to support the development of the mathematical theory of the 
constraint equations, and inspiring many others to do so. 

In this survey, we describe some of the aspects of the study of scalar curvature 
and the constraint equations that have inspired us. Asymptotically flat solutions 
of the constraints are discussed at some length, including some classic proofs of 
fundamental results, and some relatively recent work on the geometry of such initial 
data sets. We then present the conformal method for constructing solutions to the 
Einstein constraint equations based on the specification of certain components {free 
data) of the initial data set, and indicate connections to the Yamabe problem. In the 
final section, we present some results on gluing solutions of the constraint equations 
together, in such a way that the resulting manifold has initial data that resembles 
the initial data on the original pieces, and we emphasize how both conformal and 
non-conformal techniques have proven fruitful, both separately and together, in 
constructing interesting initial data sets. 

Though we include some sketches of proofs and a large collection of references, 
we do not attempt a complete treatment of any of these topics, nor do we endeavor 
to make a complete bibliography. We do hope, however, that we have included 
enough details so that readers will get a sense of some aspects of the field that have 
inspired us. 

1.0. 1. Acknowledgments. We would like to thank Adrian Butscher and Lan-Hsuan 
Huang for their comments and suggestions on an earlier draft. We also thank Hu- 
bert Bray for his support and encouragement in the preparation of this manuscript. 

1.1. Remarks on notation. Wc employ the Einstein summation convention to 
sum over repeated upper and lower indices. A semi-colon preceding an index will 
denote a covariant derivative, whereas a comma will denote a partial derivative. 
Let dvg and dfig be the volume and hypersurface measure, respectively, induced by 
a metric g, and let v or Vg be the unit normal to an oriented hypersurface. In case 
9 = gs the Euclidean metric, we let the respective measures be dve and d/Zg, 
and let the normal be Ve- We use the convention that the Laplacian is the trace 
of the Hessian. We often use subscripts to denote metric dependence on operators 
or norms, whereas such subscripts may be omitted for the Euclidean metric, so 
for example, A is the Euclidean Laplacian. Finally, we note that R{g) will be the 
scalar curvature of a metric g, whereas we will use "i?" on occasion to denote a 
large radius. 

2. The Constraint Equations 

The field equation of general relativity is the Einstein equation (with cosmolog- 
ical constant A), 

(1) Ric(g)-ii?(g)g + Ag = 8^r, 

where T is the energy-momentum tensor, and the field g is a Lorentzian metric. We 
say that a space-time {Ai , g) obeys the vacuum Einstein equation if it obeys the 
Einstein equation with T — 0- Given a manifold Ai, the left side of Equation ([ij, 
when expressed in terms of coordinates, forms a system of second order partial 
differential operators acting on the metric components g^,y which is quasi-linear, 
i.e. linear in the highest (second) derivatives of the metric. Indeed the system is 
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linear in the second derivatives of g^i, and quadratic in the first derivatives of g/ij/, 
with coefficients which are rational functions of the components g^j/. Hence the 
vacuum Einstein equation constitutes a second order system of quasi-linear partial 
differential equations for the metric g. In this discussion the manifold A4 has been 
given, whereas from the evolutionary point of view, which we adopt here, solutions 
of the Cauchy problem for Einstein's equation yield space-times A4 diffeomorphic 
to M X M, where M is an n~dimensional manifold carrying initial data, i.e., the 
initial data determines the the space-time topology and differential structure. For a 
recent survey of Mathematical Relativity, which places the questions of initial data 
into context within the study of the Einstein equation and the Cauchy problem, we 
refer the interested reader to [39] . 

An initial data set for a vacuum space-time consists of an n-dimensional man- 
ifold M together with a Riemannian metric g and a symmetric tensor K. In the 
non- vacuum case we also have a collection of non-gravitational fields which we col- 
lectively label J' (usually these are sections of a bundle over M). The vacuum 
constraint equations express the vanishing of the normal (to M) components of 
the Einstein equation, and can be expressed in terms of the initial data using the 
Gauss-Codazzi equations. If g is the metric induced on a space-like hypersurface 
M in a (time-oriented) Lorentzian manifold (A^,g), we let R^jke be the curvature 
tensor of g. If we let Kij be the second fundamental form of M in A^, and let 
Ol^jki be the space-time curvature tensor, the Gauss-Codazzi equations provide the 
following relationships: 



Here n is the unit time-like normal to the hypersurface, and the Latin indices refer 
to component directions tangent to the hypersurface M. 

Contractions of ([2])- ([3]), along with the Einstein equation, allow one to express 
the Einstein constraint equations in the following form, where we have allowed for 
the presence of non-gravitational fields: 



where R{g) is the scalar curvature of the metric g, J = —T{n, •) is the momentum 
density of the non-gravitational fields, p = T{n, n) is the energy density, and C {J-, g) 
denotes the set of additional constraints that might come from the non-gravitational 
part of the theory. The first of these equations is known as the momentum constraint 
and is a vector equation on M . The second, a scalar equation, is referred to as the 
scalar, or Hamiltonian, constraint, while the last are collectively labeled the non- 
gravitational constraints. As an example, for the Einstein-Maxwell theory in 3+1 
dimensions, the non-gravitational fields consist of the electric and magnetic vector 
fields E and B. In this case we have p — jdi^l^ -I- J — {E x B)g, and we 

have the non-gravitational constraints (for vanishing charge density) C{E,B,g) ~ 
{div g E , div g B) = 0. From the initial data point of view, we shall regard (|4|)-(l5]) 
as the definitions of the quantities J and p. Equations (l4l)-([6]) are what we shall 
henceforth call the Einstein constraint equations, or simply the constraint equations. 



(2) 
(3) 



(4) 
(5) 
(6) 



divgK - d{tigK) 
R{g)^2K-\K\l + {irgKf 
C(^,5) 



SttJ , 
16np 
0, 
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In order for the Einstein equation to have physical and geometric relevance, one 
needs to either prescribe T, or at least to impose some condition on it, such as 
an energy condition. Of particular significance to us here is the dominant energy 
condition, which requires that 

(7) TfiiyX^Y" > for all future-directed causal vector fields X and Y. 
At the level of the initial data, the dominant energy condition becomes 

(8) P>\J\9, 

where p and J are defined in (ll])-©. One checks that the condition ([7]) holds 
on {Ai , g) if and only if ([8]) holds relative to each spacelike hypersurface in A4 . 
Note that for maximal {tTgK = 0) initial data, with a non-negative cosmological 
constant, this implies that the scalar curvature is non-negative, R{g) > 0. This in 
part accounts for why we will encounter this geometric condition so often in what 
follows. 

Equations Q-dSl) fo rm an underdetermined system of partial differential equa- 
tions. In the classical vacuum setting of n = 3 dimensions, these are locally four 
equations for the twelve unknowns given by the components of the symmetric ten- 
sors g and K. A particular case of interest is the time- symmetric [K = 0) case, for 
which the vacuum constraint equations become simply the requirement that metric 
g has constant scalar curvature, R{g) = 2A. In this paper we will focus primarily 
on the vacuum case with a zero cosmological constant. 

3. A Tour of Asymptotically Flat Solutions 

Although many of the results below extend to higher dimensions, for clarity we 
will focus in this section on the three-dimensional case (n = 3) of the constraint 
equations. 

One often models isolated gravitational systems by space-times satisfying the 
Einstein equation which admit exterior regions where the metric approaches the 
Minkowski metric at some rate. In this section we will explore initial data used to 
model isolated systems. At the level of the initial data, an isolated system may be 
modeled by a space which approaches the Euclidean metric near infinity. To be more 
precise, suppose M is a three-manifold with a compact subset C <Z M for which 

k 

M\C = U Ejn, where the Em are pairwise disjoint, and each diffeomorphic to R'^\ 

m— 1 

< 1}. Then we say that (g, K) is asymptotically flat with decay rate q provided 
each Em admits coordinates for which we have — 5ij){x)\ = 0(|x|^l"l^'') 

and \d^Kij{x)\ = 0(|a;|-l''l-i-«), for |a| < ^-|- 1 and \(3\ < £, where £ G Z+ wiU be 
chosen depending on the problem at hand. For simplicity of presentation, we can 
take q = 1, though q > ^ generally gives sufficient decay for our purposes here, cf. 
[5]; see also the remark on the decay rate following the Positive Mass Theorem in 
Section 13.31 below. 

A family of explicit solutions of the vacuum Einstein equations which has proved 
important not only for physics but also for geometry is the family of Schwarzschild 
space-times. These space-times are characterized by rotational symmetry. There 
are coordinates (t, x) in which the Schwarzschild space-time metric takes the form 
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where 5b is the Euchdean metric. The parameter m is cahed the mass of the 
space-time. The space-hke shce t = is asymptoticahy flat and conformaUy flat 
with vanishing scalar curvature, and the metric gs{x) = (l + j^^) qe extends to 
a complete metric on the set \ {0}, with two asymptotically flat ends. We will 
refer to this Riemannian metric as the Schwarzschild metric below. The two-sphere 
a;| = Y inside this slice is totally geodesic, and the three- manifold has a reflection 
symmetry across it. This minimal sphere is called the horizon of the time symmetric 
slice. In the Schwarzschild black-hole space-time itself, this horizon is the central 
leaf of the three-dimensional null hypersurface comprising the actual event horizon. 

Another important family of solutions of the vacuum Einstein equations is given 
by the Kerr space-times, which are expressed in Kerr-Schild form as 

+ a^z^ V r / 

where ~ —dt^ + dx^ + dy^ + dz"^ is the Minkowski metric, and 

1 z 

9ndx^ = dx'^ — -35 7 [r{xdx + ydy) + a{xdy — ydx)] — —dz , 

r^ + a r 

where f is defined implicitly as the solution of the equation 

~4 ~2,'2|2|2 2\ 22 n 

r — r (x + y + z ~ a ) — a z — . 

This space-time models a rotating isolated system, with intrinsic angular momen- 
tum ma along the z-axis. 

In the past few decades, there has been a tremendous amount discovered about 
asymptotically flat solutions to the constraint equations. One of the most striking 
connections was used by Schoen in his resolution of the Yamabe problem. In- 
deed, if {M,g) is a closed three-manifold with positive scalar curvature R{g) > 0, 
then for any p £ M there is a positive Green's function Gp for the conformal 
Laplacian: (— Ag -I- ^R{g))G = 6p. Let G = AnGp. In an appropriate coordinate 
system y about p (conformal normal coordinates, cf. [SI])) we have an expan- 
sion G{y) = \y\~^ -I- ^ -|- 0(|y|), where ^ is a constant. Now consider the blowup 
(M = M\{p},g = G^g). By the formula for the scalar curvature under a conformal 
change of metric, we see that R{g) = 0. A direct calculation shows that in these 
coordinates, gij{y) — [(1 + 4A|?/|)(5ij + 0(|?/p)] , so that changing coordinates 
by the Kelvin transform x — we have 

g.,{x) = \y\'g.M = (l + ^) % + Oi\x\-'). 

Thus we see g is asymptotically fiat with vanishing scalar curvature. The constant 
A is related to the energy of the Ricci-flat space-time which has {M,g) as a totally 
geodesic Cauchy hypersurface, as we review below. 

In the next few subsections, we will discuss normal forms near infinity for asymp- 
totically fiat solutions of the constraints, in both the time-symmetric and general 
cases. We will also survey the relation between geometry and the mass and the 
center of mass of asymptotically flat solutions. 

3.1. Weighted spaces. Suppose {M,g) is an asymptotically flat three- manifold. 
It has proven effective to use weighted Sobolev and Holder spaces W_f{M,g) and 



6 



JUSTIN CORVINO AND DANIEL POLLACK 



C^p{M,g), respectively, to capture asymptotics of functions and tensors near in- 
finity. Let cr > 1 be a smooth function wliicli equals near infinity in an asymp- 
totically flat chart, and let 7 be a multi- index. A weighted L^-norm (p > 1) is then 
given by 



= / {\u\a'-Ya-^ dvg 
Jm 



More generally, weighted Sobolev spaces W_f{M,g) are given by 



h\<k h\<k 



M 



1/p 



We note that the weighting convention is not universal (cf. [33], [H]), and that 
we have chosen to follow [5], which uses the re-scaled measure dvg in three 
dimensions (and, of course, cr~" dvg in n dimensions). 

We now recall the basic definitions of weighted Holder spaces on {M,g). Let a 
be as above, let a{x,y) — min(cr(a;), tT(y)), and let d{x,y) be the geodesic distance 
with respect to g. We define the weighted Holder seminorm on functions / by 



One can naturally extend this to tensor fields, using parallel transport. The 
corresponding weighted Holder spaces C^'^ (M, g) are given by those functions 
/ g C'^'"(M) so that the following norm is finite: 



... xeM 

7|<fc 



The Sobolev embedding allows one to turn integral estimates into pointwise 
estimates. Indeed in dimension three, we have that if fc — - > a > 0, then ||u||„o,c. < 

C||M||^yfc,p, and in fact \u{x)\ — o(|a;|~*) as \x\ 00, cf. [5^. (The same holds in 
n dimensions in case k — ^-^ > a > 0.) For large fc, then, one can get supremum 
estimates on derivatives of u as well. 

We note that asymptotically flat metrics may also be defined in terms of these 
spaces, by requiring coordinates in an end for which the components (gij — Sij) 
belong to a weighted space, with respect to a Euclidean metric on the end. 

One reason it is nice to work with the weighted spaces is that the Laplace op- 
erator Ag is Fredholm (except for a discrete set of weights), and in an appropriate 
weight range, it is an isomorphism. We have the following weighted elliptic esti- 
mates and regularity for the Laplace operator (for simplicity, in three dimensions) 

Proposition 3.1. Let r, 77 > 0, a S (0, 1), g > n, 1 < p < q, 5 eR\Z, and k>2. 

Suppose (M, g) is asymptotically flat in the sense that in appropriate coordinates 
in any end, g is uniformly equivalent to gE, and {gij — 5.ij) G W^^^'"^ . Then 
Ag : W_f WZ^_f is Fredholm; if S e (0,1), this map is an isomorphism, and 
there is a C > so that ||w||Ti^fc,p < CII A„i(;||,,,fc-2,p 



Suppose {gij — Sij) G ' . Then for 6 G (0,1), Ag : Cjg C' 



-<5~2 



ts an 



isomorphism, and there is a C > so that \\w\\^k,a < C\\Agw\\^k-2.a . Moreover, 
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if w G C^jj and AgW e C'l^^'" D , then w G C'l'" , and there is a C > so that 



We remark that the weights that are avoided in the above result are precisely 
those that correspond to rates of growth of Euclidean harmonic functions on ]R'^\{0} 
(so that for dimension n > 3, the exceptional weights are given by (—6) £ {m G Z : 
m <2 — n, or m > 0}). 

3.2. Harmonically fiat asymptotics. Consider the time-symmetric case K — 
of the vacuum constraints, which reduces to the vanishing of the scalar curvature. 
A particularly simple form for the geometry on an asymptotic end is for g to be 
conformally flat with vanishing scalar curvature, which Bray has named harmoni- 
cally flat. Suppose that g — u'^gs, in suitable coordinates gij{x) — u'^{x)Sij, then 
the scalar curvature R{g) — —8u~^Au is non-negative if any only if Au < 0, and 
the scalar curvature vanishes if and only if u is harmonic. Any harmonic function 
u which tends to 1 as |a;| — > oo admits, by way of composing u — 1 with the Kelvin 
transform x i— >■ (inversion in the unit sphere), an expansion in spherical 

harmonics, which we write 



The number m in this expansion is called the ADM mass [H [S], and a simple 
calculation shows 



This limit exists and defines the ADM mass m{g) for the more general asymptoti- 
cally flat metrics we defined above when R{g) £ L^{M), for instance in case (g,K) 
solves the vacuum constraints. The mass plays an important role in the geometry of 
positive scalar curvature, as well as in the physics of isolated gravitational systems. 

Schoen and Yau [107j showed that metrics which are harmonically flat on the 
ends are dense in the space of asymptotically flat metrics with non-negative scalar 
curvature. 

Proposition 3.2. Suppose {M,g) is asymptotically flat with R{g) > 0. For any 
e > 0, there is a metric g within e of g (in a weighted norm), with R{g) > 0, and 
which is harmonically flat near infinity in each end E, with \m{g) — Tn{(j)\ < e. 
There is also a metric g with R{g) = which is harmonically flat at each end, with 
m{g) < m{g) + e. 

Sketch of proof. Let < < 1 be a smooth cutoff function so that '4>{t) = 1 for 
t < 1 and tp{t) = for t > 2. Fix tq, and on each end choose asymptotically 
flat coordinates defined for |a;| > rg. For 6 > tq, let tpgix) = ip{\x\9~^); ipg 
extends smoothly from the ends to all of M. Now consider the metric geix) = 
ipe(x)g{x) -|- (1 — 'ipg{x))gE{x). This metric is identical to the Euclidean metric for 
|a:| > 29, but the scalar curvature may fail to be non-negative for 9 <\x\ < 29. We 
now use a conformal deformation to change back to non-negative scalar curvature. 
Indeed, we have R{u^ge) = —u~^{8AggU — R{gg)u). We want to impose that 
R{u'^ge) = tpeR{g)] 4'eR{g) is compactly supported, and \^gR{g)\ < C9^^^'^ — 
C9^^ (in case q = 1). Let u—l + v and set E:{v) — — u^^(8Ag^M — R{ge)u). Then 



(9) 



+ 



m = lim 

167r R^+ao 




{\x\ = R} 
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E: is smooth in v, for v in C^Jg{M), or W^f{M) {p > 3/2, 5 G (5, 1)). If we let DS 
be the hnearization about v = {u = 1), we get DS{w) = — SAggW — AR{gg)w. 
Now, R{ge) is compactly supported and (for sufficiently large 9) of small norm, 
so that DS is a small perturbation of an invertible operator (the Laplacian), and 
thus it is an isomorphism. Moreover, the inverse norm is bounded uniformly in 
6] therefore by the inverse function theorem, for 9 large enough, we can solve 
E:{v) = ipgR{g) > 0. The solution v will be close to zero, so that u > and g = u^g 
will be close to 5, as desired, and the masses wiU be close, cf. |107[ 1100] . 

Suppose instead we wish to deform to zero scalar curvature. For this we study 
the linear operator (Ag — ^R{g)). We now recall the Schoen-Yau argument that 
this operator is invertible on weighted spaces, say W'^g as above, for g asymp- 
totically flat with R{g) > 0, in fact more generally in case R{g) has suitably 
small negative part R[g)^ . The operator (Ag — ^R{g)) is Fredholm of index zero, 
since it is a compact perturbation of the Laplacian. For any w in the kernel, we 
have by integration by parts (using the decay of w) and the Holder inequality, 
\\^M\h{dv,) ^ 4Ri9)~\\L^/^dv,)\\M\l2'(^d^^y where 2* = 6 is the Sobolev conju- 
gate exponent to 2, in dimension three. For R{g)~ small in L^^^, we see that w 
must be zero by the Sobolev inequality: ||w||i2* (^^^^-^ < C\\Vw\\]^2(^j^.^^-f, and therefore 

Ag — ^R{g) : W^f — >■ W^f_g is an isomorphism. 

Now we consider the second claim in the proposition, where we wish to push the 
scalar curvature to zero and arrange harmonic asymptotics. We may employ the 
result of the preceding paragraph to arrange zero scalar curvature. Indeed, we let v 
be the solution of AgV — ^R{g)v = ^R{g), which is equivalent to AgU—^R{g)u — 0, 
with u = 1 + V tending to 1 at infinity. Note that the constant functions — 
and w+ = 1 are sub- and super-solutions, respectively. By the maximum principle 
and Harnack inequality, cf. |58| , < m < 1. Since u also admits an expansion 
in any end given by u{x) = 1 + j|j + 0{\x\^'^), cf. |S], we see A < 0. Thus the 
metric u'^g is asymptotically flat with vanishing scalar curvature, and moreover, 
m{u^g) = rn{g) + 2A < m{g). 

Thus to prove the second claim and arrange harmonic asymptotics with the 
given mass bound, we may assume without loss of generality that R{g) = 0. The 
argument in the flrst paragraph of the proof can now be applied to flnd u > so 
that R{u'^gg) = 'ipgR{g) — 0, for large enough 9. The result follows. □ 

Bray [TB] observed that the Schoen-Yau approximation could be modified to 
produce ends that are precisely Schwarzschild, preserving non-negative scalar cur- 
vature. We now recall this argument. Let : R — R be a smooth, non-negative 
function with support spt((/!)) = [—1,1], which is constant near the origin. Let 
f{x) = be the associated smooth, rotationally symmetric bump function 

supported in the unit ball, and by scaling we may assume / <p(x) dve = 1. Let the 

R3 

family (pt{x) = e^'^ip{^) for e J, be the corresponding approximate identity; note 
that iff: has unit integral, and is supported on {|a;| < e}. We can use ip^ to mollify 
functions by convolution: ((^g * w)(x) = J ipe{y)w{x — y) dve- 

J/GR3 

We now state and prove Bray's proposition. The closeness of the metrics g and 
g can be measured in a norm, or, as Bray states it, as an e-quasi-isometry, i.e. for 
aU nonzero v e TM, 4^ e {e~'e'). 
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Proposition 3.3. Suppose {E,g) is an asymptotically flat end with R{g) > 0. 
For any e > 0, there is a metric g with R{g) > 0, which is e-close to g, which 
is isometric to a (Riemannian) Schwarzschild metric near infinity in E, and for 
which \'m{g) — m{g)\ < e. 

Proof. By applying Proposition I3.2[ we may modify tlie metric and choose co- 
ordinates on E and an tq > 0, so tlrat for |a;| > ro, tlic metric Iras the form 
9ij{x) — u{x)'^6ij, with Au = 0, and u{x) = 1 + + 0(|a;|^^). Now for any 
R> ro and S > 0, consider the harmonic function v{x) = Ci + with Ci and C2 
chosen so that Ci + %■ — max u(x) + S, and Ci + ^ — min u(x) — 6. By the 
choice of Ci and C2, we have that the function w defined by 



is continuous. Furthermore, note that by the expansion of u, we have that for any 
rj > 0, if we take R sufficiently large and S — ^ with a sufficiently small, then 

|Ci — 1| < ?7 and IC2 — ™2^^ | < rj. Since u and v are harmonic, and the minimum of 
harmonic functions is (weakly) superharmonic, we have that w is superharmonic. If 
we convolve w with a spherically symmetric mollifier ips supported in {y : \y\ < 6} 
as above, we produce a smooth, superharmonic function w — ^ps which satisfies, 
by the mean value property, wix) = u(x) for \x\ < R — S, and w{x) = v{x) 
for \x\ > 2R + S. Thus if we let (jij = w^Sij on E, then g agrees with g for 
To < |a;| < i? — S, and g is precisely Schwarzschild on |x| > 2R + (5, with mass 



3.3. Positive mass and topological obstructions to positive scalar curva- 
ture. In the late seventies, Schoen and Yau had the tremendous insight relating 
the proof of the positivity of the mass of an asymptotically fiat metric of non- 
negative scalar curvature to topological obstructions to positive scalar curvature 
(PSC) on closed manifolds |103| . Lohkamp [S3] later observed an interesting way 
to further link these ideas, which we describe below. We first recall the celebrated 
fundamental group obstruction to positive scalar curvature from Schoen- Yau [10^. 

Theorem 3.4. Suppose M is a closed, orientable three-manifold. Suppose further- 
more that either 7ri(M) contains a finitely-generated non-cyclic abelian subgroup, 
or that 7ri(M) contains a subgroup abstractly isomorphic to the fundamental group 
of a closed Riemann surface of positive genus. Then M admits no metric with pos- 
itive scalar curvature, and any metric on M having non-negative scalar curvature 
is flat. 

We now discuss some elements of the proof. Recall that if S is a closed orientable 
smooth minimal hypersurface with unit normal v and second fundamental form 
j4 = (hij), then the second variation of area for variation V — (piy is 




\x\ < R 
R<\x\< 2R 
\x\ > 2R 



m{g) = 2C1C2 « m{g). 



□ 



(10) 



^Cipdfig:^ / {\Vj:ip\^ ~ {\\A\\^ +mCg{l^,iy))ip^))dHg 



E E 



where Cip = A^ip + (||A|p + RiCg(;^, v))ip is called the Jacobi operator on S. 
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One of the key components of the proof of Theorem [33] is the following beautiful 

observation of Schoen and Yau, using the stability inequality — J ipC{(p) d^g > 0. 

s 

Proposition 3.5. Let {M,g) be a closed, oriented Riemannian three-manifold with 
positive scalar curvature. Then (Af , g) admits no stable minimal immersion S S-^ M 
from a closed Riemann surface E with positive genus. 

Proof. Suppose S is a closed, oriented surface, with a stable minimal immersion 
to M. Choose a local orthonormal frame {Ei, E2, E^} adapted to S, with Ei, E2 
tangential, and E3 — v, the oriented normal to E. Let hij ~ g(\7 EiEj , E3) denote 
the second fundamental form (z, j = 1, 2) of E. Let Kij denote the curvature (in M) 
of the Ei-Ej two-plane, and let denote the Gauss curvature of E, which is half 
the scalar curvature of E; we also have RiCg{E3, E3) = K13 + K23, and the scalar 
curvature R{g) — 2{Ki2 + if 13 + K2'i). The stability inequality for the variation 
field V — ipv, for S C^{M), is written (in the local frame) 

2 

j {Ki3 + K23 + J2 ^Ij)^'' '^^'9 < j dfig. 

Now one may combine minimality, hn + ^,22 = 0, along with the Gauss equation 
= K12 + ft.11ft.22 — fti2i to yield K^. = K12 — ^J2^ii- Putting this into the 



stability inequality, with ip = I, we obtain 



2 ^ '"ij- 



J [-R{g) - + 2 E ^^ff ^ 



i.j=i 

The Gauss-Bonnet Theorem J Kj^ dfig = 27rx(E) then implies 



E 



2 



(11) < i / {Rig) + ^ ft2^.) d^,g < 27rx(E). 

Thus x(S) > Oi so E has genus zero. □ 

To prove Theorem l3.41 then, one seeks to prove the existence of a stable immersed 
surface of positive genus, given the condition on the fundamental group. For exam- 
ple, given an abelian subgroup of rank two in 7ri(M), one gets a continuous map 
of a torus into M which maps 7ri(T^) onto this subgroup as follows: consider a 
torus as a rectangle with opposite sides suitably identified, and map these opposite 
sides to generators of the rank-two subgroup; the boundary of this rectangle maps 
to a nuU-homotopic curve in M, and so the continuous map extends to the interior 
of the rectangle, and hence to the torus. (A similar procedure works for higher 
genus g by representing the surface as a suitable quotient of a ig-gon.) Amongst 
all maps inducing the same action on the fundamental groups (a conjugation may 
be invoked to keep track of the base point), one finds an energy- minimizing, hence 
harmonic, map. The energy is defined with respect to a surface metric, and is con- 
formally invariant. By varying across conformal classes of (or the higher genus 
surface), one finds a map of least energy whose action on 7ri(M) is the same as the 
original. The energy-minimizer can be shown to be a stable minimal immersion. 
We can now apply the preceding proposition. See |105| for details. Note that by 
PT|) . we see that one may allow R{g) = 0, for example in case {T^,g) were the flat 
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torus. By the scalar curvature deformation results of Fischer- Marsden [5T], g must 
be flat, else the scalar curvature may be made positive. 

Now we recall an observation due to Lohkamp. We give a somewhat simpler 
proof than in [83 , using Bray's Proposition 13.31 



Proposition 3.6. Suppose the end {E,g) is harmonically flat with negative mass. 
Then there is a metric on E which has non-negative scalar curvature which is not 
identically zero, which agrees with g near dE and which is flat outside a compact 
set. 

Proof. Without loss in generality, we may assume by applying Bray's result that 
the metric g is Schwarzschild near infinity in E: in appropriate asymptotically flat 
coordinates, we write gij{x) = (1 + ^7|)^<^ij, with m < 0. For R > — y sufficiently 
large so that E contains a neighborhood of = R in coordinates, we consider the 
positive continuous function U on E given by 



[/(x)=min(l + ;^,l + ^)=r ' 



\x\ < R 

2|x|'^ ■ 2RJ 1^1 + ^ >i? 

Since U is the minimum of two harmonic functions, it is superharmonic. We use 
a spherically symmetric moUifier as above to mollify U: let U = (tpe *U). U \s 
actually well-defined on E for e small enough, and it is smooth and positive. Indeed 
for < R—e, U {x) = u{x), by the mean value property of harmonic functions, and 
for \x\ > R + e, U{x) = 1 + ^ is a positive constant. The mollification preserves 
superharmonicity, and thus, as it is clear that U is not harmonic near \x\ = R, there 
is a region where AU < 0. Using the chosen coordinates, the metric gij = If^Sij has 
non-negative scalar curvature which is not identically zero, and it is a fiat metric 
outside |a;| > i? + e. □ 

A basic version of the Positive Mass Theorem follows as a corollary of this result, 
using the Schoen-Yau topological obstruction to positive scalar curvature. 

Theorem 3.7 (Riemannian Positive Mass Theorem). Suppose {M,g) is an asymp- 
totically flat three-manifold with non-negative scalar curvature R{g) > 0. Then the 
ADM mass of any end is non-negative. 

Proof. From Proposition 13.21 we may assume that R{g) = on M, and that g is 
harmonically flat on the ends. In fact, we now recall how we can further reduce to 
the case where there is only one asymptotically flat end, cf. [171 P- 202-205] . 

If there are several (fc > 1, say) ends of M, choose one, say E, and let m be a 
harmonic function, AgU = 0, with u{x) ^ 1 as |a;| — oo in E, and u{x) ^ as 
— > oo in the other ends. To find such a u, we fix G C°°(M) with w = 1 
near infinity in the end E, and w = near infinity in the other ends. Then 
AgW £ C^{M), and so we can solve AgV = —AgW, for w in a weighted space (so 
that V decays to zero in each end). Let u = v + w. We could also use constant 
super- and sub-solutions u+ — 1 and m_ = and use a barrier argument. By the 
maximum principle, < u < 1 on M. 

Since R{g) = 0, the metric u^g has vanishing scalar curvature, and near infinity 
in any end, it can be written u'^g = U^gE for J7 > 0, where AU = 0. In the 
ends other than E, U tends to zero at infinity, and so if we write U in spherical 
harmonics in these ends, we have U{x) = jfy + ■ ■ " • The higher spherical harmonics 
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are not everywhere-positive, so that we see that c > 0. A simple calculation (again 
using the Kelvin transform x ^ x\x\^'^) shows that (M, u'^g) can be completed to a 
smooth asymptotically flat manifold {M,g) by adding fc — 1 points, corresponding 
to compactification of all but the chosen end. Since m < 1, we have the mass m{g) 
is less than that of {E,g), as in the proof of Proposition 13.21 

If the ADM mass of {E, g) were negative, then the ADM mass of (M, g) would 
also be negative. Now we apply the preceding proposition to assert the existence of 
a metric on M with non-negative scalar curvature which is flat outside a compact 
set. We can thus consider a region C M so that dW = {x : |a;'| = ro, i = 1,2, 3} 
is a large cube in the region where the metric is flat. The metric thus descends to 
a metric with non-negative (not identically zero) scalar curvature on the quotient 
space M obtained by identifying the opposite coordinate faces in pairs. M can be 
expressed as a connected-sum of a closed three-manifold with the torus T'^, and in 
particular there is a copy of Z Z in 7ri(M). We clearly have a contradiction to 
the Schoen-Yau obstruction to positive scalar curvature. □ 

The argument above yields the fact that a metric on M.^ which has non-negative 
scalar curvature and is Euclidean outside a compact set must in fact be globally 
flat. Schoen and Yau in fact prove a strong rigidity statement: if the mass of any 
end vanishes, then {M,g) is isometric to (R'^, ijb). From this we see that the decay 
rate to the Euclidean metric of an asymptotically flat metric with non-negative 
scalar curvature is constrained by the mass: for instance, any such metric which 
has the expansion in asymptotically flat coordinates gij{x) — Sij + 0{j^) with 
q > 1 must in fact be flat. Finally, we note that the negative mass Schwarzschild 
metrics, given by gsix) = (1+ 2pT)*-9^' with m < 0, have vanishing scalar curvature 
and an asymptotically flat end. However, they become singular at \x\ ~ > 0; 
which is at finite geodesic distance, and such metrics are incomplete. 

The minimal hypersurface proof of the Positive Mass Theorem given by Schoen 
and Yau jl04j does not use Theorem 13.41 directly, but the proofs share some of 
the same ideas. A calculation in asymptotically flat coordinates for a metric 
g.jix) = (l+^)%-fO(|a;|-2) yields the Christoffel symbol rf^. = ir^6,,+Oi\x\-^). 
By tracing over i,j — 1,2, we see that if the mass were negative, the coordi- 
nate hyperplanes \x^\ = a for large enough a have mean curvature vector H with 
g{H, gfa) < 0. Hence these hyperplanes can be used as barriers for finding a 
stable minimal hypersurface asymptotic to a plane, by solving a Plateau problem 
on large cylinders with axis along the x'^-direction. The stability inequality and 
Gauss-Bonnet yield a contradiction, in a similar manner as above, cf. [lOOj . 

3.4. Global charges and Killing initial data (KIDs). Space-times which are 
sufficiently asymptotically flat possess conserved quantities which are sometimes 
called global charges, one example of which is the ADM mass we have already met. 
In the framework of Noether's Theorem, conserved quantities correspond to symme- 
tries. Generic space-times may not possess symmetries, but in asymptotically flat 
space-times, the background Minkowskian Killing vector fields provide good enough 
approximate symmetries to induce global conserved quantities, as measured by a 
Minkowskian observer (coordinate chart) at infinity. These charges and their cor- 
responding fields are as follows: mass-energy, linear momentum, center of 
mass, x^-§f — t-^', angular momentum, — x''-^. The center of mass and 
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angular momentum fit together to form a two-form J^,y, whose components give 
the charge associated to the Minkowski Killing field x^^ — x'^ [TH 137] . 

We find it convenient to use the momentum tensor tt*^ = K"^^ — {tTgK)g^^ , for 
which the vacuum constraint equations become 

(12) + i (tr^^)2 _ |^|2 ^ 

(13) divgTT = 0. 

We write the constraints map as $(5, tt) = {R{g) + 5 (tfgTr)^ — |7r|g, divgTr), so that 
the vacuum constraints correspond to solutions $((7,7r) = (0,0). 

The relationship between linearization stability of the Einstein equation, the 
kernel of the formal adjoint of the linearization and the existence of 

space-time symmetries forms an important backdrop for results that will follow. 
We refer to the works of Fischer-Mardsen 52J, Moncrief [93 , and the more recent 
works [inillS] for more details. We assume (M, g, tt) satisfies the vacuum constraint 
equations, and consider £)$ at (5,7r). Moncrief showed that space-time symmetries 
(Killing fields) correspond precisely to elements of the kernel of D(f>*. To illustrate, 
we note that if g is the Ricci-flat metric determined by the Cauchy data on M, and 
if ^ is a time-like Killing field for ^, so that M is a level set of t with unit time-like 
normal n and adapted coordinates , we can write ^ = Nn + X = Nn + g|j , 
and so g = -N'^dt^ + gij{dx^ + X''dt){dx^ + XHt). Thus {N,X) forms the lapse 
and shift for the stationary metric g. It can be shown [93j that D^*{N ,X) = 0. 
An element {N,X) in the kernel of £)$* at [g^Ti) is called a KID [Killing Initial 
Data). Conversely, a non-trivial KID {N,X) for {M,g,TT) can be used to generate 
a (not necessarily time-like) Killing field in the Einstein evolution of the Cauchy 
data iMlIIS]- 

At the Minkowski data {gE,0), D^*{N,X) = {-{AN)gE + 'ilessN,-^LxgE), 
where Lxgs is the Lie derivative. The kernel K of is thus the direct sum of 
Kq = span{l ^} together with the space of Killing fields oi gs, which of 

course is spanned by the generators of rotations and translations. 

As it turns out, asymptotic flatness alone is not enough to guarantee a well- 
defined angular momentum, cf. |151 166] . Regge and Teitelboim j98i proposed as- 
ymptotic conditions sufficient to guarantee a well-defined angular momentum and 
center of mass for solutions to the constraints, namely that in suitable asymptoti- 
cally flat coordinates the following estimates also hold (with the decay rate g = 1 
for simplicity): 
(14) 



>S{g^,{^)-g^A-^))\=0{\x\-^-^-'), \d^^[K,,{x) + K,,{-x))\=0{\x\-^^'^-'). 

These conditions impose approximate parity symmetry on the data {g, K) in an 
asymptotically flat end. The Kerr space-times admit coordinates (such as Kerr- 
Schild) for which the constant time slices (and corresponding boosted slices) satisfy 
(|T4]) . The space of initial data satisfying ([M]) is known to be dense in the space of 
vacuum asymptotically flat data |47j , as we discuss in the next section. 

Using such a coordinate system, we can compute the energy and linear and 
angular momenta using flux integrals at infinity. We let Y(j) be the Euclidean 
rotational Killing fields, e.g. F(3) = x^-^ — x^-^, we let be the Euclidean 
outward normal, and recall that the Einstein convention is in force. Furthermore, 
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let the integral over Soo denote the limit of integrals over = i?} as i? — >■ +00. 
The flux integrals for the energy and momenta are given by 



1 

T, 

S. 



So. 



(15) Pi = j [Kij - K^igij)vldHe = ^ j TTijVldUe 

Soo 

{Kjk - K^e9jk)Y^)i^idtJ,e = ^ j 3kY(;)4dtie 

(16) mc'' = ^ f [Y^x'' {gij,i - guj) - ^ {gieS'^^vl - giiU^)\ dfie- 

a i i 

Joo 

In the last term, one may replace gab in the integrand by {gab — Sab)- Taken together, 
these give a set of ten Poincare charges associated to the end. We remark that, from 
the physical point of view, the quantity m, which corresponds to the approximate 
symmetry generated by might more properly be termed the energy E instead 
of the mass, as the charges {E,P^) are the components of the energy-momentum 

four-vector. Prom the point of the view of the space-time, this vector is Lorentz- 

3 

covariant, and the square of the rest mass is £^ — ^ (P*)^. We stick to the notation 

i=l 

above for consistency with the Riemannian case. 

It is instructive to view these charge integrals as arising from integrating the 
constraint functions against elements of the cokernel of the linearized constraint 
operator. By linearizing at the Minkowski data (^sjO), we have + h.Tr) = 

D^^hyir) + Q{h,Tr), where a simple expansion (with q = 1) yields an estimate 
of the quadratic error term Q{h,Tr) = 0(\x\^^). In Euclidean coordinates at the 

Minkowski data, D^{h, tt) = ( '^Zi^ij^j^^a.n)- S '^ij-j ) • Thus for any KID [N, X), 

^ i,3 j 

i.e. for any vector and scalar pair {N,X) which satisfies D^*{N,X) = (0,0), we 
have as a consequence of integration by parts that 

(17) j {N,X)-<^>{gE + h,n)dve=B{R)-B{Ro)+ J iN,X) ■ Q{h,n)dve 

A(Ro,R) A{Ro,R) 

where A{Ro, R) = {Rq < \x\ < R}, and 

B{R) = / (mjX' + N V(/iij,i - hiij) - y2{N,ihij - Njhu))ui rf/Xe- 

J{\x\=R} ^ i i ' 

By letting X be a Euclidean Killing vector field, or letting A'' be a constant or 
a coordinate function , we can easily relate B{K) to one of the above surface 
integrals defining the ADM energy-momenta. 

3.4.1. Comparison to Newtonian theory. We now compute the mass and center 
of mass in a harmonically flat end, and compare to classical Newtonian theory. 
We expand the harmonic conformal factor u{x) = 1 + 1^ + + 0(|a;|^'^), with 
P-x = (3^x^-\-l3'^x'^-\-p^x^, and we note the error term picks up an extra |x|~^-decay 
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with each successive derivative. Using this in the flux integrals for mass and center, 
we obtain 



1 
16^ 



{\x\=R} * 



{\x\ = R} 

Thus we see m^2A and mc^ — 2 Ac'' = 2/3''. 

We note that if we translate asymptotically flat coordinates y'"' — + , then 
the center of mass integral becomes ra{cJ^ + a^) as expected. In the harmonically 
flat case, we note that the expansion transforms under translation as u{y — a) — 
1 + + ^ + 0(|y|~^), where ji^ = fi^ + Aa^ . So the dcflnition of corresponds 
to the translation which makes the |a::|^^-terms in the expansion vanish: u[y + c) = 
1 + f + 0(|y|-3). 

In the Newtonian setting, we can study gravitational systems in terms of the 
mass density p, with decay conditions on p corresponding to the system being 
isolated; for example, we might take p to be compactly supported. This gives rise 
to a gravitational potential (p which satisfles A0 = Anp (we have taken Newton's 
constant G — 1). In case p is compactly supported, cf) is harmonic near infinity; we 
take 4'{x) — >■ as |x| — > oo, and thus can be expanded as a series in spherical 
harmonics, which yields (j){x) — — — + 0{\x\~^). Note that the total mass 

J]jj3 p dve and the moments /jja x^ p dve defining the center of mass of the system 
can be written as fiux integrals in terms of the potential, analogous to those above: 

p dVf. = — / A(/) dVf. — — lim ( — ^ dp,. — m 

f 1 /" 1 /" 86 x^ 

/ x^p dve^ — / x''A(j) dVe = — lim / (x''- (f>)dpe = 13''. 

7r3 47r7js,3 47r ii^+oo7{|^|^^j ^ dr R 

We see then that the center of mass c'^ can be defined when m 7^ by mc'^ = (3^ . 

The density p generates the potential 4>. Analogously we can consider perturba- 
tions of the Euclidean metric, say j = gE + h, where h is sufficiently small and for 
simplicity has compact support. We generate from this a potential function u > 
tending to 1 at infinity, for which R{u^^) — 0, i.e. A^u = ^R{'-^)u. Near infinity 
u is harmonic, so that for large \x\^ u{x) = 1 + + + 0(|a;|~^), and we have 
(compare to the Newtonian formulas above) 

m — lim / -— dpe = / A-yU dv^ — / R(j)u dv^ 

27rfl.^+cx> J dr ^ 2tt J ^ ^ 167r ^ ^ 

{\x\ = R} R3 

/3'= = I y'^A^udv^^-^ [ y''R{j)udv^; 

47r Jr3 6ZTT Jr3 

in the preceding integral, we have taken global harmonic coordinates for the 
metric 7 (so that A-yy* = 0, which simplifies Green's formula above), with |x — ?/| = 
0(|x|-i/2), valid for h small [3 Theorem 3.1]. 
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3.5. Harmonic asymptotics. In this section we discuss an extension of harmon- 
ically flat asymptotics to the setting of the full constraint equations, in which we 
study both the metric g and second fundamental form K of the space-like slice. 
These special harmonic asymptotic conditions were proposed by Schoen in con- 
junction with a proof of the Positive Mass Theorem, cf. [102j , but have already 
proven useful in other contexts [47l |64] • 

Let CgX = Lxg — divg{X)g, where {Lxg)ij — Xi-j + Xj-i is the Lie derivative 
operator. Note that tr g{CgX) = —divg{X), in contrast to the trace-free conformal 
Killing operator T?g{X) = Lxg — ^divg{X)g that is used in the conformal method 
for solving the constraints (see Section 2]). The harmonic asymptotic conditions 
require that outside a compact set there exist a positive function u which tends to 
1 at infinity, and a vector field X which tends to at infinity, so that g = u^gs and 
TT = u^Cgj^X, and so that in suitable asymptotically fiat coordinates g.^ = u^Sij 
and TTij — v?{Xij + Xj^i — X^Sij), where X^ = Xi. For such data, the vacuum 
constraint equations become (with C — Lgj^ , A = Ag^ ) 

(18) 8Am-u(- |/:X|2 + i(tr(£X))2) =0 

(19) AX, + 4u-iuj(£X)^' - 2u-^u^,tx{LX) = 0. 

From these equations one can immediately see the advantage of using Lg: its di- 
vergence has leading order term which is precisely the Laplacian of the components 
of X. This fact is used in the proof of the following density theorem [4 7) . 

Theorem 3.8. Let S G (^il) and p > |, and suppose (A/, g, tt) is a vacuum 
initial data set with (gij — 6ij,Trij) G W^f x W^f_g. Given any e > 0, there is a 
vacuum initial data set (g, 7f) with harmonic asymptotics and within e of {g, tt) in 
the W_f X W_f_ _g norm, and so that the mass and linear momentum of {g, w) are 
within e of those of [g, vr). 

The proof is obscured, in comparison with that of the Schoen- Yau harmonic 
asymptotic approximation sketched above, by the fact that the linearization at 
a general asymptotically fiat solution of the constraints is somewhat complicated, 
as is that of the constraint map for harmonic asymptotics (the left-hand side of 
(fT8l) - (fT9|) above). The method of proof uses the surjectivity of D$ in appropriate 
spaces, and has already been employed in other works, cf. [TPl IM]. 

An important feature of these asymptotic conditions is that ((T4)) holds, and 
moreover, the total energy-momentum and angular momentum are directly encoded 
in the asymptotics of (u, X), and thus these conserved quantities directly affect the 
asymptotic geometry. In fact, one solves for {u ~ 1,X) in a weighted space, and a 
standard expansion for the solution of Poisson's equation, cf., e.g., [S], yields the 
following expansions near infinity: 

(20) u{x)^l + ^ + Oi\x\-^), X'{x)^^ + 0{\x\-^). 

\x\ \x\ 

These leading coefficients are, up to constant factors, the ADM mass/energy and 
the linear momentum of the end. Indeed, it is easily shown as above that A = m/2. 
A quick computation yields 



SCALAR CURVATURE AND THE EINSTEIN CONSTRAINT EQUATIONS 



17 



Thus on a large coordmate sphere \x\ ^ R with outward Euclidean normal i^e = ify , 
we have niji/i — T^ij^ — — ^p + Odxl^"^). Thus we see the linear momentum vector 

Though in general we cannot expand (u,X) further in spherical harmonics, we 
can in fact expand the odd parts of u and X to one higher order. From the 
flux integrals for angular momentum and center of mass, we see that the terms 
at the next order of the expansion that will contribute to the limiting values of 
the integrals are precisely the odd parts. From the constraint equations ([T8| -(fT9 |) 
for {u,X), using ([20]), we see Au°'^'^ = 0{\x\-^) and AX°'^'^ = 0{\x\-^), so that 
u°dd(-^) ^ ^_|.0(|a;|-3) g^^^ (X*)°dd(2,) ^ ^m^ + o(|a;|-3)_ We have seen how the 

vector /3 is related to the center of mass, and similarly, the angular momentum can 
be written in terms of dj^), e.g. J3 = 5(^^(1) ~ '^(2))- ^^^'^ remark the following: 
if we shift the coordinates to y = x — c, let X{y) — X{y + c) and expand, we get 

{X'Y'^'^iy) = + 0(|y|~3), where J^-j = - d'B^ = d\^ + 2c^P\ Thus we 

see that the angular momentum transforms as Ji ^ Ji — {c x P)i, as expected. 

3.6. Geometry of the center of mass. A natural question is whether the notion 
of center of mass can be described in a geometrically interesting and intrinsic way. 
This may be particularly intriguing in the vacuum case, where the only energy 
comes from the gravitational field (metric). In this section, we discuss approaches 
to answer this question, and the relation between them. 

3.6.1. Center of mass and conformal symmetries. As we noted above, the mass 
and center of mass are related to asymptotic symmetries of the space. Based on a 
suggestion of R. Schoen, L.-H. Huang [64] employed a computation of the mass and 
center using approximate conformal symmetries, as we now recall. The connection 
is motivated in part by the following generalized Pohozaev identity from |101| . 

Proposition 3.9. Let {M,g) be a compact Riemannian n-manifold with smooth 
boundary dM . Suppose X is a conformal Killing field on M . Then we have the 
following identity, where v is the outward unit normal to dM : 

(21) I X{R{g))dv, = -^j (Ric(3) ~ ^9)[X, ,.)df,, . 

M dM 

To illustrate, we note than in the case of harmonically flat asymptotics, the 
mass and center can be associated to exact conformal Killing fields in the asymp- 
totic region. Since dilations and inversions are conformal isometrics of the flat 
metric, we can represent dilation or translation of infinity by conjugating dilations 
or translations with the Kelvin transform to obtain the following conformal isome- 
trics near infinity, for any a > 0, and any y e R"^ (computed in the harmonically 
flat coordinates x, for which gij = u^Sij): 
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By taking the derivative of the first map with respect to a at a = 1, we see that 
the infinitesimal generator of this family of isometrics is the conformal Killing field 
^X(o) = —x^-^. Differentiating the second map at y = yields the conformal 
Killing fields X^t) = kPgfi - ix^x'^, for ^ = 1, 2, 3. 
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We compute the integrals of RiCg(X(£), ly) over a large sphere {|a;| = R} in the 
asymptotically flat region. We take the y-normal pointing toward infinity, and 
let = u^i^ be the Euclidean normal, so that 

f - 

RiCg(X(o), i^)d^g = -2 / —x^x^Rjkdiie = lQT:m+0{R^^), 

J \x\ 

{\x\=B} {\x\=R} 

RiCg{X(^(>),iy)dfig = 327r/?^ + 0{R-^) = 167rmc^ + 0{R-^) , 

{\x\=R} 

where we used Rjk = {Ric{g))jk = —u~^(2uujk + Quju^k — 2^jfe|Vu|g) (see [iS] 
for details) . If g is harmonically flat near inflnity, then by applying (PT|) to a large 
annular region, we get an asymptotic conservation statement for the mass integral 
across the boundaries of the region. We note that in the center of mass integrals, 
the leading-order terms from the expansion are not absolutely integrable, but they 
have the correct parity so that their surface integrals vanish. 
More generally, the ADM mass can be expressed as 

(Ric(g) - ii?(g)5)(X(o),z^) dfj.g, 

{\x\=R} 

and following [64], a center can be deflned, when m ^ and ([14]) holds, by 
(22) ^^f = Y^,H-^ / iRH9)-lRi9)9m.„.)d,g. 

{\x\=R} 

3.6.2. The center of mass foliation. In an asymptotically flat chart for {M,g), the 
coordinate spheres near infinity are approximate solutions to a constant mean cur- 
vature (CMC) equation. One can perturb these to exact solutions either by the 
mean curvature flow due to Huisken and Yau 70 , or by an implicit function theo- 
rem method due to Ye [114j , cf. [92] . The methods produce the same foliation near 
infinity in the case of positive mass, and this foliation provides a geometric notion 
of the center of mass. 

In the Huisken- Yau approach, one studies solutions F"' : §^ x / — M of the flow 
■§fF'^ = {h — H)v^ with initial condition Fq :— -f^lsa^io} the standard embedding 
of the sphere of radius a centered at the origin. Here u is the outward unit normal, 
H is the inward mean curvature of the surface F^(EP) = F'^(S^ x {t}) =: , 
and h is the integral average of the mean curvature: Jj^^,y{h. — H)dfig = 0. This 
normalized mean curvature flow improves the isoperimetric ratio (area decreases 
while the enclosed volume is fixed). For a > ao (ctq large), the solution exists 
for afl times t > 0, the surfaces -^/^(S^) converge to surfaces ^^(S^) := , and 
the limiting configuration {Af^, cr > (Tq} forms a foliation by stable constant mean 
curvature spheres. In fact, we have the following result from |70j . 

Theorem 3.10. Consider an asymptotically flat three-manifold {M,g), and sup- 
pose in an exterior region the metric can be written gij{x) — (l + ^^)^(5ij-|-0(|x|~^) 
with m > 0. There is a ao > 0, positive constants Ci and C2, and a vector 
Chy € so that for each a > ao, the following are true. The initial value problem 
has a unique smooth solution for all times t > 0. The surfaces = F^{§'^) con- 
verge exponentially fast to a smooth stable hypersurface Af^, with constant mean 



m = lim 

167r R^+00 
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curvature Ha. The radial coordinate r ~ \x\ satisfies |r — crj < C'l on M"^ , and 
\Ha — |- + ^1 < C2fT^^ ■ The hypersurfaces M'^ have a joint center of mass vector 
Chy, in the following sense (where dj^i^ is the pullback of the Euclidean surface 
measure on M" ): 

J xdfie J FS, dfia 
Chy = lim ^^"p , = lim ^ ^ . 

tr-H-oo J cl^e ct->-+oo j 

Huisken and Yau argue that the foliation is asymptotically round on approach 
to infinity, and they obtain a certain uniqueness result for the foliation. Ye likewise 
obtains asymptotic estimates of the geometry of the leaves of the foliation, and 
a uniqueness theorem. We state following result, which follows directly from the 
approach of Ye |114) : in particular, notice how the CMC foliation is produced by 
solving for a shift r and a normal perturbation of large coordinate spheres. 

Theorem 3.11. Consider an asymptotically flat three-manifold with an exterior 
region which satisfies the conditions in Theorem \3.1(A Let < a < 1 . There is a 
(To > large enough, and a constant C > Q so that for cr > ao, M'^ is the image 
of the embedding $p : §^ 9 w i— )■ p{t{p) + uj + ip{p,uj)iy{uj)) , where ^ — ^ — H^, 

■)llc2.°(s2) < ^ and t{p) < ^. If gp is the induced metric on Af^, then as 
p — > +0O, p^^^*{gp) converges in C^'" to the unit round metric on S^. 

We note here that the mass is assumed to be positive, but the constraint equa- 
tions are not imposed, and in particular, no local energy condition is assumed. 
Actually, Ye's result only requires that the mass parameter be nonzero. Positivity 
of m, however, implies that the leaves of the foliation are stable for the isoperimet- 
ric problem. It is an interesting question to what extent the leaves of the foliation 
solve a more global optimization problem. The relation of the mass to the geometry 
of the three-manifold via the isoperimetric problem arises in the work of Bray [16] 
(cf . |18| , [461 [1] ) on the Penrose Inequality, and more recently on an isoperimetric 
approach to the definition of mass of isolated systems due to Huisken [68] . 

We now have three notions of center: c from the flux integrals defining global 
charges ([16]) . Ci from ([22]) . and Chy from the geometric foliation by CMC spheres 
near infinity. We now discuss the relation between these. 

3.6.3. Equivalence of centers. We first state the Corvino-Wu result, which estab- 
lishes, in the special case where the metric is sufficiently harmonically flat near 
infinity, that the center of mass from the ADM formulation agrees with the geo- 
metric center Chy, a question raised by X. Zhang [116] (and others). 

Theorem 3.12. Consider an asymptotically flat three-manifold with an exterior 
region that admits an asymptotically flat chart in which gij{x) = (1 -|- + 

^-^^-t^^j^±^)4 5,.j+0{\x\-^), with m > 0. Then m this chart, C^y = ^ c*^. 

Sketch of the proof. By translation, it suffices to prove that in case gij{x) = (1 -I- 
2^)^(5ij -f 0(|a;|~^) with m > 0, then Chy = in this chart. The metric gij{x) — 
Ifxl^^^ij +^^(^1^^) agrees with Schwarzschild in this chart to one higher order 
in \x\~^ than is generally considered in [70] . Hence, the radial coordinate spheres 
Fq{S'^) = {\x\ = a} in this chart are better approximate solutions of the CMC 
equation; we have arranged this by design by centering them appropriately. It 
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suffices, then, to show that the flow moves them sufficiently little so that the center 
of mass determined by the foliation is the zero vector. To prove Theorem l3.12| then, 
we need to estimate carefully how much varies from Fq , and to keep track of the 
variation in the induced surface measure on . Both of these quantities evolve by 
an equation involving {H — h), as dtF"' = {h — H)v^ and dt{dfj.g) = H{H — h)dfig. 
Basic estimates from |W yield |Vi/p < Ca~^, which can be integrated to prove 
\H — h\ < Ca~^. This is what one expects for a finite center of mass. The goal, 
then, is to show \H — h\ has better decay when we choose coordinates for which 
gij = (1 + 2]7|)^'^*i + C){\x\^^). Since we are starting the flow well-centered, we 
expect the spheres to drift very little, in fact on average not at all. Indeed, we have 
the following proposition from [48]. 

Proposition 3.13. There is a C > so that for all a > ao, and for all t > 0, 

(23) max \H - h\ < Ca~'/^e-^^'''/^ . 



We can conclude the theorem by inserting this into the definition of Chy, and 
by using the flow and the initial conditions to estimate F^. □ 

In two recent papers, L.-H. Huang has definitively unified and extended the above 
notions of center of mass to the general setting of Regge-Teitelboim asymptotics 
([T4l) . which we denote "AF-RT." In [64], she shows that c = C/, and in the case of 
strong enough asymptotics for which the Huisken-Yau proof yields a foliation, that 
c = C/ = Chy ■ In [65.^ she proceeds obtain a foliation near infinity by CMC spheres 
under AF-RT asymptotics, essentially unique for m ^ 0, and which is stable for 
m > 0. The center of mass of the foliation again agrees with C/. We state the main 
existence theorem from [65] . 

Theorem 3.14. Assume (M, g, K) is an AF-RT end with decay rate g e (5: 1] '^'^'^ 
center of mass c. If m ^ 0, then for R sufficiently large, there exist surfaces S/j of 
constant mean curvature H^.^ ~ Ti.''^ 0{R^^^'^) . S/j is the graph of a function ipQ 
over the coordinate sphere Sr{c) — {x : \x — c\ — R}, i.e. Tifj = {a; -I- ipo{x)i'g : x £ 

Sb{c)}. Ifij^{x)^Mc + R^), then Uollc-.^ + R\\ii'o)°'''\\c--- < coR^^' ■ Thus c 
is the geometric center of the family S^. If m > the surfaces T,ji form a strictly 
stable foliation. 

We do not go into the proof, which involves very beautiful applications of elliptic 
theory, but we do note a connection to the above asymptotics discussion. Huang 
establishes and uses an extension of the density theorem for harmonic asymptotics. 
In [47], it is shown that solutions with these asymptotics are dense in the space of 
asymptotically flat solutions, with the energy-momentum four-vector approaching 
that of the original data. In Huang argues that if the data is AF-RT to start, 
then the odd part of g and the even part of K can also be made arbitrarily small, 
and thus so can the change in all global charges under the approximation, cf. [67] . 
We note that Huang also obtains interesting uniqueness results which one may 
compare with [TQl |96] . 

We will return to the issue of the asymptotics in Section 15.11 
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4. The Conformal Method 

The most successful approach so far for systematically studying the existence 
and uniqueness of solutions to (jH)-® is through the conformal method of Lich- 
nerowicz [82", Choquct-Bruhat and York '33'. The idea is to introduce a set of 
unconstrained conformal data, which are freely chosen, and find {g, K) by solving 
a determined system of partial differential equations. In this section, we treat the 
general n-dimensional case, as the theory is the same in all dimensions ti > 3. In 
the vacuum case with vanishing cosmological constant 33 , the free conformal data 
consist of an n-dimensional manifold M, a Riemannian metric g on M (with Levi- 
Civita connection D), a trace-free symmetric tensor a, and the mean curvature 
function r. The initial data (g, K) defined as 

(24) g = 0«5 

(25) K = (jj-^ia + VW) + ^(j)'^g , 

where q = and is a positive function, will then solve (HJ-lIS]) if and only if i 
and the vector field W solve the equations 

(26) divij(T>W + a)^ !!—l,j)i+^DT , 

n 

(27) - ^—R{~g)q^ + -^Aa + VWljcb-^-' -r^c^^+i = , 

q[n — 1) q[n — 1) ^ qn 

where V is the conformal Killing operator: 

(28) VWab = DaWb + DbWa " -gabDcW . 

n 

Vector fields W annihilated by V are conformal Killing vector fields, and are char- 
acterized by the fact that they generate (perhaps local) conformal diffeomorphisms 
of {M,g = (f>'^g)- The semi-linear scalar equation (|27p is often referred to as the 
Lichnerowicz equation. 

Equations ([26|) - ([27l) form a determined system of equations for the scalar- vector 
pair ((^, VF). The operator divg(I? •) is a linear, formally self-adjoint, elliptic op- 
erator on vector fields. What makes the study of the system (|26 l) - (f27)) difficult in 
general is the nonlinear coupling between the two equations. 

The explicit choice of is motivated by the two identities 

(29) R{g) - -<p-''-'{q{n - l)Ag<j> - R{g)<j>) , 
where g — (fi^'^g, and 

(30) b'^ir^Bab) = <p-''~^D-Bab 

which holds for any trace- free tensor B. Equation ([^^ is the well known identity 
relating the scalar curvatures of two conformally related metrics, and we note that 
q = is the unique exponent that does not lead to supplementary terms 
in (EId. 

In the space-time evolution (A^"+^, g) of the initial data set (M, g, K), the func- 
tion r = trg/C is the mean curvature of the hypersurface M G Ai. The assumption 
that the mean curvature function t is constant on M significantly simplifies the 
analysis of the vacuum constraint equations, because it decouples equations 
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and ([27l) . One can then attempt to solve ([SG]) for W, and then solve the Lichnerow- 
icz equatfon ([27]). 



4.1. Existence via the conformal method: CMC initial data. Existence and 
uniqueness of solutions of this problem for constant mean curvature (CMC) data 
has been studied extensively. For compact manifolds this was exhaustively analysed 
by Isenberg [71] , building upon a large amount of previous work |5U [M] 11151 [55] ; 
the proof was simplified by Maxwell in [84] . If we let 



(31) yi[g])^ _ inf 



/ec~(Af)j^o 11/11 



2 



where 2* = 2n/{n — 2), denote the Yamabe invariant of the conformal class [g] of 
metrics determined by g (see [81]), the result reads as follows: 

Theorem 4.1 (|71|). Consider a smooth conformal initial data set {g,a,T) on a 
compact manifold M , with constant t. Then there always exists a solution W of 
(|26p . Setting a ~ VW + a, a positive solution cj) of the Lichnerowicz equation exists 
if and only if one of the following conditions holds 

1. y{[g])>0,a^Q. 

2. 3;([5]) = 0,(T^0,r^0. 

3. :y([5]) <0,T^0. 

4. y{[g])^0,a = {),T = 0. 

Moreover, we have uniqueness in every case except case 4 (where one may homo- 
thetically scale the metric to choose the volume arbitrarily). 

Beyond a simple existence result, we see that Theorem 14 . 1 1 does more. It actually 
solves the Conformal Parametrization Problem for CMC initial data on a compact 
manifold. In particular, given a smooth, compact Riemannian manifold {M,g), 
it provides a full parametrization of the set of vacuum solutions (g, K) to the 
Einstein constraint equations such that g lies in the conformal class [g\ of g. This 
parametrization is given by the choice of a constant r and a trace-free symmetric 
tensor a such that one of the four conditions of Theorem 14.11 are satisfied (note 
that a vanishes if and only if a vanishes). The extent to which the Conformal 
Parametrization Problem may be solved for non-CMC initial data is unclear, and 
is currently actively being explored, as we discuss below. 



Remarks on the proof of Theorem It is worth commenting on the role of the 
Yamabe invariant in this result. Indeed the Lichnerowicz equation bears a great 
deal of similarity to the Yamabe equation specifying the existence of a conformal 
metric h — (fflg of constant salar curvature R{h) — C 

(32) A 0_^^i?(g)0+^^C</.'+i^O. 

q(n — 1) q{n — 1) 

There are two key features which distinguish the two equations: one is the existence 
of the term involving <f>~''~^ in the Lichnerowicz equation and the other is the sign of 
the term involving 0'+^. These distinctions are crucial. Indeed the solution of the 
Yamabe problem [99] [81] was one of Richard Schoen's significant contributions to 
geometric analysis and a very important part of late twentieth century mathematics. 
The original proof of Theorem 14.11 used the solution of the Yamabe problem in full. 
Maxwell showed in [84] that the proof in fact relies only on the existence of a metric 
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of within a given conformal class whose scalar curvature has the same sign as that 
of the Yamabe invariant 3^([5])- This far simpler fact is expressed in the following 
well known and very useful result. 

Proposition 4.2. Let Lg = — Ag + ^(„"'_i) -^(ff) ^6 the conformal Laplacian, and 
let fig be its first eigenvalue. Then the following are equivalent. 

1- > (respectively — or < 0). 

2. fig > (respectively = or <0). 

3. There exists a positive function u on M such that the scalar curvature of 
u'^g satisfies R{u'^g) > (respectively — or < 0) everywhere on M. 

Maxwell established this result in a low regularity setting in [53] in order to 
carry out a program for analysing the conformal method for solving the constraint 
equations for metrics of low differentiability (see also [13] )• This was motivated 
in part by work on the evolution problem for rough initial data [80\ \79\ \78\ IllOj . 
Judicious applications of this result, together with the maximum principle, allow 
one to prove Theorem 14 . 1 1 using the method of sub and super-solutions. □ 

The conformal method easily extends to the CMC constraint equations for many 
non- vacuum initial data sets, e.g. the Einstein-Maxwell system ^Tlj, where one 
obtains results very similar to those of Theorem 14. II (see also [72]) ■ However, other 
important examples, such as the Einstein-scalar field system [28j |30l [29l |6T] are not 
as well understood. 

4.2. Existence via the conformal method: near-CMC and far-from-CMC 
initial data. Conformal data which is close to being CMC (e.g. as measured 
by a smallness assumption on |VT|/|r|) are usually referred to as "near-CMC." 
Classes of near-CMC conformal data solutions have been constructed [75| , [27| , [2l,l76 ] . 
In particular, with an appropriate definition of near-CMC, it is know that on a 
Riemannian manifold with no conformal Killing fields, the conditions of Theorem 
14.11 yield a unique solution of the constraint equations provided that the mean 
curvature is nowhere vanishing. There is at least one example of a non-existence 
theorem [77] for a class of near-CMC conformal data. However, due to the nonlinear 
coupling in the system (P5)) - (P7)) . the question of existence for unrestricted choices 
of the mean curvature r appears to be significantly more difficult, and until recently 
all results assumed strong restrictions on the gradient of r. 

Alan Rendall has demonstrated that on x endowed with the product metric, 
the data consisting of a = and r equal to an (arbitrary) odd function on §^ has no 
solution which shares the symmetry of the conformal data. Therefore either there 
are no solutions, or there are more than one solution. 

The first general existence result in the "far-from-CMC" context is due to Hoist, 
Nagy, and Tsogtgerel [62] [63]. In the Yamabe positive ([V([g]) > 0) case, they 
construct solutions with freely specified mean curvature; however they assume both 
the presence of (sufficiently weak) matter fields and that a is pointwise sufficiently 
small (depending on r) and not identically zero. One may view the smallness 
assumption on a as standing in for the near-CMC hypothesis in previous work. 
In [87], Maxwell provides a sufficient condition, with no restrictions on the mean 
curvature, for the conformal method to generate solutions to the vacuum constraint 
equations on compact manifolds. As an application. Maxwell demonstrates the 
existence of a large class of solutions to the vacuum constraint equations with 
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freely specified mean curvature (again assuming that a is pointwise sufficiently 
small and not identically zero). In the results of Hoist, Nagy, and Tsogtgerel, as 
well as those of Maxwell, the methods employed do not allow one to assert the 
uniqueness of the solutions found. Nonetheless, these results together represent a 
significant advance in our understanding of how the conformal method may be used 
to generate solutions of the vacuum constraint equations. However the existence 
question for generic classes of large conformal data remains wide open. 

In order to explore the "large data" regime more systematically. Maxwell has 
recently studied a model problem in the Yamabe null (3^([g]) = 0) class [551. By 
considering a three-parameter family of model conformal data that allow for simul- 
taneous violations of both the near-CMC and small-tj conditions. Maxwell was able 
to identify a number of new phenomena. First, he confirmed that for this class of 
data there is also a small-a result. He was also able to assert the non-existence 
of solutions when the data violated both the near-CMC and small-a- conditions. 
Moreover, he was able to demonstrate that the small-a solutions are not unique. 
These results indicate that the landscape for studying the constraint equations via 
the conformal method will be quite interesting when one moves away from small 
data as represented by either the near-CMC or small-o- conditions. It seems likely 
that in order to effectively analyse this situation, a new approach to studying the 
constraint equations may be needed. 

In [111 the reader will find a presentation of alternative approaches to construct- 
ing solutions of the constraints, covering work done up to 2003. 

4.3. The constraint equations on asymptotically flat manifolds. There are 
a large number of well-established results concerning the existence of CMC and 
near-CMC solutions of the Einstein constraint equations on asymptotically fiat 
manifolds [22l [HH |33l [Ml EH ESI |86] using the conformal method. In 1977, Cantor 
introduced a quantity, analogous to the Yamabe invariant, which was sufficient to 
prove the existence of a positive solution to the Lichnerowicz equation relative to 
a given set of asymptotically flat CMC conformal data ^22j. Since in the asymp- 
totically flat setting CMC means that r = 0, the relevant Lichnerowicz equation is 
simply 

(33) A, ^ - + ^1-1^^"'^^ - « ■ 

Maxwell |85j uncovered an error in Cantor's definition of the invariant [22| and 
provided the correct definition. For an asymptotically fiat manifold {M,g), we 
define 



(34) yAF{[g])^ inf „,,,2 

(compare with (|31l) ). The precise result is then the following. 



Theorem 4.3. Suppose that {M,g) is asymptotically flat of class Wj^ for k > n/2 

and T G (0, n — 2) and that a G W^^J^ is a transverse-traceless tensor. Then there 
exists a positive conformal factor (j) satisfying hS!^) if and only if yAF{[9\) > 0. 
Moreover, if a solution exists, then it is unique. 

The key step to establishing this is to have an appropriate analog of Proposition 
14.21 This is given by the following result. 
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Proposition 4.4. Suppose that (A/, g) is asymptotically flat of class W'l'^ for k > 
n/2 and r G (0,ri — 2). Then the following are equivalent. 

1. There exists a conformal factor 4> > such that 1 — (j) £ W_^'^(Af) and such 
that (j}'^g is scalar flat. 

2. yAFM) > 0. 

3. For each rj G [0, 1], P,, = — Ag + V q^n-i) ^is) i^*^ isomorphism acting on 

WtriM). 

5. Gluing Constructions 

Gluing constructions, by wliicli known solutions of a geometric partial differential 
equation are combined to produce new solutions, are now ubiquitous in geometric 
analysis. One of the earliest gluing results concerns scalar curvature. In 1979, 
Schoen and Yau [106] showed that, within the category of compact n-dimensional 
manifolds, the property of admitting a metric with strictly positive scalar curvature 
is preserved under surgeries of codimension fc > 3. In particular, this includes tak- 
ing connected sums, which can be viewed as codimension n surgery. (Gromov and 
Lawson [591 160j independently established these results, introducing other impor- 
tant techniques into the study of manifolds with positive scalar curvature.) Over 
the past ten years, there have been many applications of gluing constructions to 
general relativity. All of these take place at the level of the constraint equations 
and have implications for space-times by considering the evolution of the initial 
data sets constructed by gluing. 

5.1. Asymptotic gluing. It was an interesting open question for many years to 
what extent the asymptotic expansion at infinity determines the interior behavior 
of a solution of the (vacuum) constraint equations, cf. |108| p. 371], [8]. The sim- 
plest such question was answered in the resolution of the Positive Mass Theorem: 
if near infinity the solution is precisely Euclidean, then the solution is globally flat. 
A natural next question to pose, then, is whether an asymptotically flat metric 
of vanishing scalar curvature, which near infinity agrees precisely with a standard 
time-symmetric slice of a Schwarzschild solution, must be globally Schwarzschild. 
The question is to some extent about unique continuation under the constraint 
of vanishing scalar curvature, inspired possibly in part due to the success of the 
conformal method discussed above, which turns the constraints, which are unde- 
iermmed-elliptic, into a determined elliptic system. We remark that it was known 
that unique continuation does not hold with respect to the interior, as Bartnik 
[9] constructed non-trivial scalar-flat metrics on M.^ with a fiat sub-domain. We 
also note that, as we recalled above in Proposition 13.31 Bray showed that there 
are plenty of metrics which are Schwarzschild near infinity and have non-negative 
scalar curvature, and in particular. Cutler and Wald [49 showed that there exist 
solutions of the Einstein-Maxwell constraints on R"^ which near infinity agree with 
a Schwarzschild metric. We do remark, however, that the Schwarzschild metric 
satisfies a rigidity condition, with respect to the Penrose inequality [171 169) . 

In the late nineties, Corvino and Schoen resolved the question in a very strong 
form, as we discuss below. The main idea which echoes in their work is that there 
is a lot of freedom in the initial data, and one ought to look outside the space 
of conformal deformations. The idea that took hold was a localized version of 
the Fischer-Marsden results (SU (52] , the possibility for which was in part inspired 
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from the work of Lohkamp [83 . It was shown that not only do there exist many 
initial data sets for the vacuum constraints which agree with Schwarzschild (time- 
symmetric) or Kerr near infinity, but furthermore such solutions are dense in an 
appropriate topology. 

We now state the time-symmetric result from |44| . 

Theorem 5.1. Let {M,g) be an asymptotically flat three-manifold with zero scalar 
curvature. Let E d M be any asymptotically flat end, and let Erg be an exterior 
region in E corresponding to {x : \x\ > tq} in asymptotically flat coordinates. Let 
k be a non-negative integer. Then for any e > 0, there is an R > and a (smooth) 
metric g with zero scalar curvature and \\g — 5||c"=(E) < ^ (t^^ norm is taken with 
respect to the Euclidean metric in the asymptotically flat coordinate chart), so that 
g is equal to g on M \ Ej^, and g is identical to an asymptotically flat end of a 
standard Schwarzschild slice on E2R. The analogous statement holds for n > 3. 

The general case of the constraints was addressed in [33, cf. [37] • We define a 
family of solutions on the exterior of a fixed ball and smoothly parametrized on an 
open set O C M^'' to be admissible if, with reference to a fixed coordinate chart 
near infinity, the family satisfies ([T4| locally uniformly, and the map Q : O ^ 
which associates to each member of the family its energy-momenta (m, P, J, mc) 
is a homeomorphism onto an open subset. We note that slices in Kerr form an 
admissible family, parametrized by the total mass to, the angular momentum pa- 
rameter a, and an element in the Poincare group to represent Euclidean motions 
of the asymptotic coordinate system as well as boosts, cf. [37] . 

Theorem 5.2. Let (AI,g,TT) be any asymptotically flat solution of the vacuum 
constraints. Given any e > 0, there is a solution (g, vf) within e of (<?, tt) (in a 
weighted norm) and whose ADM energy-momentum {E,P) is within e of that of 
{g,Tr), so that near infinity, (g, tt) agrees with a member of an admissible asymptotic 
model family, for example a boosted space-like slice in Kerr. 

We now discuss the main ideas of the proof. In a given asymptotically flat end, 
we pick an asymptotically fiat coordinate chart, and at a large coordinate radius 
\x\ — R, we patch our given data to that of the model using a cutoff function in 
the annulus Afi from |xl = i? to |x| = 2R. Since the data is approaching the 
fiat data, the gluing produces an approximate solution of the constraints. At this 
point we seek to perturb this to an exact solution. This is often done using the 
conformal method; however, the conformal factor induces a global change, albeit a 
small one far away from the gluing region, and we seek to keep the data unchanged 
outside the annulus Ar. The key that enables one to do this is the underdetermined 
nature of the constraints, i.e. the overdetermined-ellipticity of the adjoint of the 
linearized constraint operator. In fact, the construction relies heavily on linear 
elliptic estimates for this adjoint, as we now discuss. 

5.1.1. Localized scalar curvature deformation. Motivated by the preceding discus- 
sion, we study the following localized scalar curvature deformation problem: Given 
a smooth domain D, C {M,g) and a compactly contained sub-domain ^Iq, for / 
supported in and sufficiently small, find a tensor h supported in O so that 
R{g + h) = R{g) -\- f , with h small depending on /. A natural way to approach 
this, in the spirit of Fischer-Marsden [51 , is to study the linearization of the scalar 
curvature operator. 



SCALAR CURVATURE AND THE EINSTEIN CONSTRAINT EQUATIONS 27 

Let Lg be the linearized scalar curvature operator: Lg{h) — ^|j_q-R(.9 + th). 
Then we have Lg{h) — —Ag{tTgh) + divg(divgft,) — h ■ Kic{g), and thus L*u = 

— {Agu)g+}iesSg u—uRic{g). By taking the trace of this equation, we see trg{L*u) — 

— {n—l)AgU—uR{g), so that AgU, and hence the full Hessian, is controlled pointwise 
by L*u and u. So we get the following estimate, valid on any domain (without 
imposing boundary conditions): 

(35) \\u\\H2^n)<C{\\L;u\\L^n) + Mmin))- 

We can replace the lower-order 77^-term by the L^-norm on smooth bounded do- 
mains follows by using interpolation on the compact inclusions ^ ^ L^. 
In fact if we restrict u to lie in a subspace transverse to the kernel of L* (generically 
there is no kernel, but we will be interested in the case where g is the flat metric, 
which has a four-dimensional kernel) , then we have an estimate with no lower-order 
term at all. 

This estimate can be used to get analogous estimates in weighted spaces. Indeed, 
let p be a smooth, positive function on il, with p{x) = {d{x))^ near dfl, where 
d = d{x) is distance of x to d^l, and N will be taken to be sufhciently large (one 
may also take p{x) = e~^^'^^^^ near dfl); we take p to depend only on d and to 
be monotonic in d, leveling off to a positive constant away from dfl. From the 
above, choosing u to lie in a fixed subspace transverse to the kernel of L* (if it is 
non-trivial), we get an estimate ||u||^f2(Q^-) < C||L*u||i2(Q^), for small e > 0, where 
fie — {x G n : d{x) > e} and C is uniform in e. For example, such a uniform 
estimate holds for metrics near the flat metric gE, for functions u transverse to the 
kernel Kq = span{l, x^, x^, x^}. By multiplying the square of this estimate by p'{e) 
and integrating by parts, we get the weighed estimate ||w||_ff2(f2) < C'\\L*u\\]^2(^fiy 
There are similar weighted estimates for which are slightly harder to prove, 

but hold nonetheless. 

We illustrate the usefulness of such weighted estimates by presenting a simpler 
example which W. Qiu [97] adapted from the analysis in [44] . Here we look at 
the divergence operator, whose adjoint is —V, say at the flat metric. Suppose 
we want to solve divX = / on an open ball B. Suppose that / has compact 
support in i?, and that we want X to decay to zero at dB. Of course a necessary 
condition on / is that its integral over the ball vanish, Jg / dve = 0, i.e., / must be 
orthogonal to the kernel of the adjoint of div. We solve this problem variationally: 
let J-"(m) '■— (^|Vm[^/9 — uf) dve, where p is as above; for the problem at hand, 
we could follow [92] and also choose p to be identically 1 on the support of /, but 
this is not essential. We minimize J-{u) over all u G Hl^^{B) so that u^p and jVu^p 
are integrable (these integrals can be used to define the (squares of) the weighted 
and iJp-norms), and so that u is L^-orthogonal to C, where C > is a smooth 
bump function with compact support. Since u is taken to lie in a space transverse 
to the kernel of V, we have an estimate (as above) of the form 

\\u\\hI(b) < C||Vu||i2(B) , 

which is simply a weighted Poincare inequality, cf. [97 . Thus 

Standard functional analysis yields a minimizer u. The Euler-Lagrange equation is 
simple to compute, and indeed letting X = —pVu we see that divX ~ f + X( for 
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some A G M. By the elliptic regularity, we see that u is smooth in the interior of 
B, and moreover, by integrating the Euler-Lagrange equation against the constant 
function 1, we see that A = and so X solves the original equation (the integration 
by parts is justified by the pointwise estimate below). Now although u might not 
decay at dB, we have that X does, by the Schauder estimates. Indeed, let x G i? 
have distance d to dB, and let Bi and B2 be the balls of radius | and ^ about x; 
we can take c? < 1 so small that B2 is outside the support of /. We have the elliptic 
equation Au + ^ • Vu = —j;- We take p equal to a large power N of the distance 
to the boundary, near the boundary. Then by standard elliptic estimates, we have, 
with ip — 'ip{n, k, a) and using / = on B2, and J^^ u^dve < Cd~^ J^^ v? p dve, 

mic'-HB,) = \\p^u\\c...(B,) < Cd""-^ (ll/p-ic-i.»(B.) + 

< Cd^'^'\\uhi(B)<Cd^-^\\fh.^^^sy 

The last inequality follows from the inequality J^{u) < and the weighted estimate. 
From this, we see that we can solve for X decaying at the boundary along with as 
many derivatives as we like, provided we choose N large enough. 

We now return to the scalar curvature problem. Suppose that L* has trivial 
kernel on fi. Essentially the same analysis, using the weighed estimates and the 
functional J^{u) ~ J {^p\L*u\'^ — uf) dvg, allows use to solve Lg{pL*u) — f vari- 
n 

ationally. For instance, if / is small, then h — pL*u will be small enough so that 
g + h g + pL*gU is a metric, and we have R{g + h) — R{g) + / + Q{h). Thus 
we see that in case L* has trivial kernel, then we can choose / small and solve for 
the localized scalar curvature deformation, at least to first order. We can iterate 
linear corrections to solve the non-linear problem, i.e. we can solve the localized 
scalar curvature deformation problem as posed above. On the other hand, metrics 
like the Euclidean metric, for which L* has non-trivial kernel, are called static, 
and elements in the kernel are sometimes called static KIDs. Static metrics are 
very special; for instance they must have constant scalar curvature. This local 
deformation result cannot hold at static metrics; for example in connection with 
the Positive Mass Theorem, the flat three-torus does not admit any positive scalar 
curvature metric, and analogously, there is no compactly supported deformation of 
the Euclidean metric with non-negative scalar and non-constant scalar curvature. 
The round sphere is also static, and it was only recently that Brendle, Marques 
and Neves provided a counterexample to the Min-Oo conjecture: they construct 
metrics on n > 3, which have R{g) > n{n — 1), with strict inequality holding at 
some point, with g a unit round metric in a neighborhood of a hemisphere |19] . 

5.1.2. Localized deformations, global charges, and the proofs of Theorems \ 5.1\ and 
\5.SX To prove Theorems 15.11 and 15. 2[ we patch together two solutions of the con- 
straints in annulus A^i near infinity, producing an approximate solution. The data 
there is close to (g_BjO), so to get uniform estimates on the adjoint of the lin- 
earized constraint operator, we have to work transverse to Kq = ker L* , respec- 
tively K = ker D<^* . Indeed, using the methods discussed above along with suit- 
able projections, we can, at the linear level, solve the constraints up to the kernel 
Kq, respectively K . By iterating the linear solvability in a Newton-Picard itera- 
tion, we can produce g (respectively {g,TT)) which solves R(cj) G QKq, respectively 
$(5, tt) G C-^ (for C a suitable bump function), and which agrees with the glued 
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approximate solution except in the gluing region An] that is, it is exactly the model 
solution on £'27?, exactly the original on AI \ Er, and a perturbation of the glued 
solution in Ar. The key is that we have used the overdetermined-elliptic estimates 
to solve for a compactly supported (in An) perturbation. 

So we now have our constraint operator with values in (^K; of course we want 
this value to be zero. Naturally (as with the divergence example above), we take 
our constraint values and integrate against K to see what we get. For example, for 
the kernel element 1 G ker L*, we get (compare ([T7)l) 

R{g)dve^ / ^igij,ij - gii,jj) dve = / ^{gij^i - gii^j)vldiie 

Ar J Ar ^ j JOAr j 

« IGttAto, 

where Am is the mass of the outer Schwarzschild minus the mass of the original 
data. Analogously, integration against a linear function G yields the change 
167rA(mc). For the KiUing field X = ^fr G ker £»$*, we have 

/ (divgS-)jX* dve ~ / (divgg7r),X* dv^ = / tt^j-XV^ d^e « SttAPi. 

JAr J Ar JOAr 

Integrating against a rotation field gives the change in angular momentum. 

Thus we see that the constraint equations imply an asymptotic conservation 
law across the boundary of the annulus. The conserved quantity is exactly the 
boundary integral coming from the divergence theorem applied to the inner product 
of the constraint operator and the generator of the asymptotic symmetry, i.e. the 
element of the kernel of the adjoint of the linearized constraint operator. For our 
gluing procedure, we see that our model near infinity has to roughly preserve the 
conserved quantities coming from the kernel of the linearized constraints. It is this 
kernel, then, that dictates which part of the asymptotics has to be preserved in our 
construction, and so our model solutions have to constitute a family which exhibits 
all possible values of the asymptotic quantities in an effective way. Thus far, our 
gluing and perturbation procedure depended only on the asymptotic fiatness, and 
so we were not quite able to solve the constraints. We now pay attention to which 
one of the model family we have glued on, and show that one of these works to 
solve the constraints. 

In the time-symmetric vacuum case, we are trying to make the metric g be 
scalar-flat. Of course, what we actually have is a family of metrics, one for each 
of the Schwarzshild family (1 -I- 2\^c\ )'^9e we glue on, parametrized by the mass 
m and center of mass c (relative to the fixed asymptotic chart). Our procedure, 
then, produces a map from the parameter values to the values R{g) G C^Oi a map 
between four-dimensional Euclidean spaces. The map roughly looks like (up to 
scaling and constant factors) 

(36) (m, c) M> (m — toq: 'TIC — moco) , 

where mo and cq are the mass and center of the original metric; note we can always 
translate coordinates so that cq = 0. This is the leading term, and the error term 
goes to zero as R —> -\-oo. For large enough i?, then, the parameter map is a small 
perturbation of a local diffeomorphism which covers the origin at the initial mass 
and center, and thus by degree theory, the actual map hits zero for some (m, c) 
near (mo,co). For these parameter values, R{g) = 0. The analogous results hold 
for the general case. 
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In summary, we see analytically that the obstruction to doing local gluing comes 
from the cokernel of the linearized operator, and that this kernel exactly corresponds 
to quantities that must be (almost) conserved in the procedure. Thus to meet this 
constraint, we need a good family of models at infinity for which the conserved 
quantities are effective parameters. We also emphasize that we need to make good 
approximate solutions, which we can do here since the data tends to flat near 
infinity. 

Note that this general model for a gluing proceedure is well known and often used 
throughout geometric analysis. The key novel feature introduced here is to exploit 
the underdetermined nature of the constraint equations (via the overdetermined- 
ellipticity of the adjoint of the linearized constraint operator) to solve for the per- 
turbations with compact support. 

We make two further remarks. We saw earlier that if we deform the Euclidean 
metric with a small compactly supported tensor h, we can impose the constraints in 
the conformal class of + ft., and the asymptotic expansion of the conformal factor 
has a positive coefficient of the |a;|~^-term. In this way, we put the appropriate 
mass at infinity through the conformal factor. For the asymptotic gluing, the 
localized deformations do not quite solve the problem: we also need to put the 
appropriate mass (and other charges) at infinity. Rather than doing this via a 
conformal deformation, which would induce a global perturbation, we establish 
this deformation by choosing an member of our model family of solutions at infinity 
which captures this necessary deformation in these asymptotic parameters. Finally, 
we remark that this analysis has been extended to allow for the gluing of multiple 
asymptotic ends satisfying the vacuum constraints into a single asymptotic end, 
which may be interpreted as initial data for the gravitational A^-body problem 

[22135]. 

5.1.3. Asymptotically simple vacuum space-times. In the 1960's, Penrose [9 5) pro- 
posed a model of isolated gravitational systems (asymptotic simplicity) based on 
the conformal compactification of Minkowski space, cf. [112] . Indeed if we write 
the Minkowski metric as = ~dt^ + dr^ + r'^gs^ , where gs^ is the metric on a round 
unit two-sphere, we can then introduce advanced and retarded null coordinates 
V = t+r, u — t—r, and re-scaled time and radial coordinates T = arctanw-t-arctanw, 
R — arctanw — arctanu. In these coordinates we obtain, for i}-^ = l{l+v^){l+u^), 

T] ^ Q-^ {-dT^ + dR^ + svn^{R) 532) = n-'^{-dT'^ + g^s). 

Thus we have an embedding of Minkowski space-time into the Einstein static uni- 
verse R X which is a conformal isometry. The image of Minkowski space is a 
pre-compact subset, the boundary of which represents infinity, time-like, space-like, 
and null, of Minkowksi space-time. It is the essential features of the construction, 
which represents infinity conformally faithfully but at a finite distance, that Penrose 
tries to capture in his definition of asymptotic simplicity. 

The Penrose proposal has had enormous influence on the study of gravitational 
radiation, so one would naturally like to establish the existence of a rich class of 
such space-times, for instance through stability results which yield new examples 
through perturbation. Friedrich addressed the stability problem by rewriting the 
Einstein equation to emphasize the conformal structure, and he obtained a small 
data, semi-global stability result: for hyperboloidal data suitably close to a given 
hyperboloidal data set in Minkowski space (intersecting future null infinity), the 
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resulting solution of the initial-value problem for the Einstein vacuum equation 
admits a conformal compactification to the future [54] , cf. [55l \57\ and the more 
recent work of Anderson and Chrusciel [3]. On the other hand, work by Friedrich 
[56] and more recent work by Valiente Kroon [111] show that the structure of initial 
data near space-like infinity must indeed be special for the evolution to admit a 
smooth conformal compactification. 

The stability results of Friedrich provide a way to construct non-trivial asymp- 
totically simple space-times: control the asymptotics near spatial infinity on an 
asymptotically fiat initial data set, in such a manner that the data will evolve to 
a space-time with suitable hyperboloidal slices (to the future and the past), then 
apply the stability result to evolve from here. In fact. Cutler and Wald [49] use 
this method to produce examples of such space-times for the Einstein-Maxwell 
equations. They construct a family of time-symmetric initial data (generated by 
scaling a specially constructed magnetic field), which approaches Minkowski data 
as the mass tends to zero, such that outside a fixed ball each member of the family is 
identically an end of the standard asymptotically fiat initial slice of a Schwarzschild 
space-time. This still left open the question of whether there exist non-trivial purely 
radiative space-times, i.e. non-trivial vacuum space-times (with A = 0) which admit 
a conformal compactification in the sense of Penrose. In recent years this question 
has been resolved, and in fact, we have the following theorem. 

Theorem 5.3. There exists an infinite- dimensional family of solutions (R^, g) of 
the time- symmetric vacuum constraint equations ( with A = 0^ whose Einstein evo- 
lution is an asymptotically simple space-time, i.e. the maximal Ricci-fiat space-time 
(M^,g) with the three- geometry (R"^,;;) as a totally geodesic Cauchy surface admits 
a conformal compactification in the sense of Penrose. 

To prove this, one shows the existence of families of solutions of the time- 
symmetric vacuum constraints approaching the fiat metric which can be perturbed 
to nearby solutions that are Schwarzschild outside of fixed radius. The difficulty 
which arises is that in applying the methods of [44], the resulting center of mass 
may drift outward as the mass tends to zero, cf. ([36]) . and one must get estimates 
to control this. Chrusciel and Delay [36] first proved this statement by applying 
the methods of [44j to produce an infinite-dimensional space of parity symmetric 
solutions; under parity, the center of mass is always zero. A different approach is 
taken in [45) , where the idea is to construct a large family of time-symmetric initial 
data sets on R'^ with a harmonically flat end, tending toward the Euclidean metric, 
for which the higher-order asymptotics are suitably bounded in terms of the mass 
TO. This allows the required estimates to scale down with m tending to 0, while 
retaining an estimate on the center of mass. 

We briefly sketch the argument, details for which can be found in 45j. We 
remark that the construction not only yields examples, but establishes what may 
be interpreted as a (weak) stability result for asymptotic simplicity at the Euclidean 
metric. First recall the York decomposition [1151 [22]: a symmetric tensor h which 
decays at inflnity (in a weighted space) can be written (for an appropriate vector 
field X) as h = h'^'^ -\-'D{X)-i'^ti{h)gE, where h^^ is transverse-traceless (trace-free 
and divergence free), and V is the conformal Killing operator introduced earlier, 
and we compute at the Euclidean metric. Let /i be a compactly supported solution 
of the linearized scalar curvature constraint L(h) = at the Euclidean metric, so 
that the TT-part hP"^ from the York decomposition is non-trivial. We note that 
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a TT tensor with respect to the flat metric is in the kernel of L, and it is known 
that there is an infinite-dimensional space of compactly supported TT tensors at 
the flat metric, cf. [HISSIIS^- If we let = u^{gE + ^h), and m{e) = m{gg), then 
L{h) — implies m'(0) = 0. We can estimate the mass m{e) from below, following 
the Brill-Deser approach to the local positivity of the ADM mass [20 {cf. 32 ). 
Indeed, the second variation of the mass is given by 167rrn"(0) = 5 / \Vh'^'^\'^ dve- 

This allows us to estimate m(e) > e^, while the fact that R{gE + e^) ^ yields a 
weighted estimate of := {u^ — 1), for large — ^ ^ m{e). This 

is the key estimate to control the center of mass. 

5.2. Conformal gluing constructions. In [73], Isenberg, Mazzeo and Pollack 
developed a gluing construction (often referred to as "IMP gluing" ) for initial data 
sets satisfying certain natural non-degeneracy assumptions. The perspective taken 
there is to work within the conformal method, and thereby establish a gluing the- 
orem for solutions of the determined system of PDEs given by (f26| and ((27|l . This 
was initially done only within the setting of constant mean curvature initial data 
sets and in dimension n = 3 (the method was extended to all higher dimensions 
in [72]). The construction of [75] allowed one to combine initial data sets by taking 
a connected sum of their underlying manifolds, to add wormholes (by performing 
codimension-3 surgery on the underlying, connected, three-manifold) to a given 
initial data set, and to replace arbitrary small neighborhoods of points in an initial 
data set with asymptotically hyperbolic ends. 

In [90| . building on work he had done for constant scalar curvature metrics 
[89l [91], Mazzieri generalized the IMP gluing to the setting of the Schoen-Yau 
and Gromov-Lawson surgery result discussed earlier. Namely he showed that un- 
der certain conditions, CMC solutions of the vacuum constraint equations may be 
combined into new solutions by gluing along a common (i.e. isometrically embed- 
ded with diffeomorphic normal bundles) submanifold of codimension greater than 
or equal to 3. 

In [74] the IMP gluing construction was extended to only require that the mean 
curvature be constant in a small neighborhood of the point about which one wanted 
to perform a connected sum. This extension enabled the authors to show that one 
can replace an arbitrary small neighborhood of a generic point in any initial data 
set with an asymptotically flat end. Since it is easy to see that CMC solutions of 
the vacuum constraint equations exist on any compact manifold |113| . this leads 
to the following result which asserts that there are no topological obstructions to 
asymptotically flat solutions of the vacuum constraint equations. 

Theorem 5.4 ( 74J). Let M be any closed n-dimensional manifold, and let p S M. 
Then M \ {p} admits an asymptotically flat initial data set satisfying the vacuum 
constraint equations. 

5.3. Initial data engineering. The gluing constructions of [T^ and [T^ are per- 
formed using the determined elliptic system provided by the conformal method, 
which necessarily (due to the unique continuation property for this system) leads 
to a global deformation, small away from the gluing site, of the initial data set. 
Now, the ability of the Corvino-Schoen asymptotic gluing technique to establish 
compactly supported deformations invited the question of whether these conformal 
gluings could be localized. This was answered in the affirmative in [37] for CMC 
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initial data under the additional, generically satisfied ^14 . assumption that there 
are no KIDs (cokernel of the linearized operator, of. Section [3^ in a neighborhood 
of the gluing site. Namely, Chrusciel and Delay showed that, assuming the absence 
of KIDs, an additional perturbation can be made to localize the IMP conformal 
gluing. This left open the question of whether generic, localized gluing could be 
performed from the start, with making any of the assumptions (like CMC) which 
were present in the IMP gluing from the dependence on the conformal method. 

In HOI HI], this question was answered, and the IMP gluing was substantially 
improved upon, by combining the gluing construction of [73 together with the 
Corvino-Schoen gluing technique of [44l [37], to obtain a localized gluing construc- 
tion in which the only assumption is the absence of KIDs near points. For a given 
n-manifold M (which may or may not be connected) and two points pa G M , 
a = 1, 2, we let M denote the manifold obtained by replacing small geodesic balls 
around these points by a neck §"~^ x /. When M is connected this corresponds 
to performing codimension-n surgery on the manifold. When the points pa lie in 
different connected components of M, this corresponds to taking the connected sum 
of those components. 

Theorem 5.5 ( |40l 141] ). Let {M,g,K) be a smooth vacuum initial data set, with 
M not necessarily connected, and consider two open sets fia C M, a = 1,2, with 
compact closure and smooth boundary, such that the set of KIDs within each VL^ 
is trivial. Then for all G i^a, e > 0, and A; G N there exists a smooth vacuum 
initial data set {M ,g{e),K{e)) on the glued manifold M such that {g{e), K(e)) is 
e-close to {g,K) in a x topology away from B{pi,e) U B{p2,(). Moreover 
{g{e), K{e)) coincides with {g,K) away from VLi UJ72- 

This result is sharp in the following sense: first note that, by the positive mass 
theorem, initial data for Minkowski space-time cannot locally be glued to any- 
thing else which is non-singular and vacuum. This meshes with the fact that for 
Minkowskian initial data, we have non-trivial KIDs on any open set $7. Next, recall 
that by the results in |14], the no-KID hypothesis in Theorem [53] is generically sat- 
isfied. Thus, the result can be interpreted as the statement that for generic vacuum 
initial data sets the local gluing can be performed around arbitrarily chosen points 
Pa- In particular the collection of initial data with generic regions Via satisfying the 
hypotheses of Theorem 15.51 is not empty. 

The proof of Theorem 15.51 is a mixture of gluing techniques developed in [72] [73] 
and those of [47] [44] [37]. In fact, the proof proceeds initially via a generalization 
of the analysis in [75] to compact manifolds with boundary. In order to have CMC 
initial data near the gluing points, which the analysis based on [73] requires, one 
makes use of the work of Bartnik [6] on the plateau problem for prescribed mean 
curvature spacelike hypersurfaces in a Lorentzian manifold. 

One application of Theorem 15.51 concerns the question of the existence of CMC 
slices in space-times with compact Cauchy surfaces. In [7], Bartnik showed that 
there exist maximally extended, globally hyperbolic solutions of the Einstein equa- 
tions with dust which admit no CMC slices. Later, Eardley and Witt (unpublished) 
proposed a scheme for showing that similar vacuum solutions exist, but their argu- 
ment was incomplete. Using Theorem 15.5] we obtain: 
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Corollary 5.6 ([40[ I41j). There exist maximal globally hyperbolic vacuum space- 
times with compact Cauchy surfaces which contain no compact spacelike hypersur- 
faces with constant mean curvature. 

Compact Cauchy surfaces with constant mean curvature are useful objects, as 
the existence of one such surface gives rise to a unique fohation by such surfaces [H] , 
and hence a canonical choice of time function (often referred to as CMC time or 
York time). Foliations by CMC Cauchy surfaces have also been extensively used in 
numerical analysis to explore the nature of cosmological singularities. With this in 
mind, one natural and important question is the extent to which space-times with 
no CMC slices are common among solutions to the vacuum Einstein equations with 
a fixed spatial topology. It is, in particular, expected that the examples constructed 
in [m HO] are not isolated. 

5.3.1. Non-zero cosmological constant. Gluing constructions have also been car- 
ried out with a non-zero cosmological constant @3j |42j |38] . These constructions 
yield space-times which coincide, in the asymptotic region, with the corresponding 
black hole models. For time-symmetric slices of these space-times, the constraint 
equations reduce to the equation for constant scalar curvature R = 2A. Gluing 
constructions have been previously carried out in this context, especially in the 
case of A > 0, but in [43l |42l |38] the emphasis is on gluing with compact support, 
in the spirit of the Corvino-Schoen technique. 

The time-symmetric slices of the A > Kottler-Schwarzschild-de Sitter space- 
times provide "Delaunay" metrics (see [43] and references therein), and the main 
result of [43l |42] is the construction of large families of metrics with exactly De- 
launay ends. When A < the focus is on asymptotically hyperbolic metrics with 
constant negative scalar curvature. With hindsight, within the family of Kottler 
metrics with A G R (with A < corresponding to the Schwarzschild-anti-de Sitter 
metrics, A = the Schwarzschild metrics and A > the Schwarzschild-de Sitter 
metrics), the gluing in the A > setting is technically easiest, while that with A < 
is the most difficult. This is due to the fact that for A > 0, one deals with a single 
linearized operator with a one-dimensional kernel, which corresponds to changing 
the Delaunay parameter; in the case A = 0, the kernel is {n -\- 1) -dimensional cor- 
responding to changes in the mass and center of mass; while for A < 0, one needs 
to consider a one-parameter family of operators, each having (n -|- l)-dimensional 
kernels [38] . 
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